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SUMMARY 

A method of developing an input data file for use in the aerodynamic 
analysis of a complete airplane with the VORLAX computer program is described. 
Familiarity with VORLAX input data and panelling methods is assumed. The 
required hardware consists of an Interactive graphics terminal equipped with a 
graphics tablet. Software Includes graphics routines from the Tektronix PLOT 
10 package as well as the VORTAB program described herein. 

The user must determine the size and location of each of the major panels 
for the aircraft before using the program. Data is entered both from the ter- 
minal keyboard and the graphics tablet. The size of the resulting data file 
is dependent on the complexity of the model and can vary from ten to several 
hundred card Images. After all of the data has been entered, two programs 
READB and PLOTB, are executed which plot the configuration allowing visual 
inspection of the model. 


INTRODUCTION 

To describe an aircraft configuration using the VORLAX aerodynamic computer 
program (ref. 1 ) a set of geometric data arranged in a specific format is 
required. At times it may be necessary to analyze an aircraft configuration 
for which the required numerical data is either nonexistent or in a form that 
is incompatible for use with the VORLAX program. In most cases, however, an 
aircraft three-view, drawn to scale, is available. Since a three-view drawing 
contains most of the required geometric data, it is a logical starting point 
from which to develop a numerical model. 

Program VORTAB has been written to generate a VORLAX input data file 
from information obtained from the configuration three-view drawing. VORTAB 
permits the definition of as many major and minor panels as required to define 
the configuration within the limitations imposed by the VORLAX program. The 
types of panels that can be defined are cambered and uncambered flat panels, 
fuselage panels, and nacelle panels. Camber data cannot be determined with 
sufficient accuracy from the three-view drawing and must, therefore, be entered 



in tabular form. The aircraft configuration is assumed to be symmetrical about 
its X-Z plane. 

The procedure file PVORTAB (used to execute program VORTAB) and program 
VORTAB have been written to be used on the Langley Research Center CDC computers 
with the Network Operating System (NOS) and Tektronix Plot 10 graphics package. 
Langley Users may find the procedure file and computer program on access files 
under user number 214737C. 

PROGRAM DESCRIPTION AND SAMPLE EXECUTION 

A description and sample execution of each section of the program will 
now be presented, illustrated by copies of the terminal display during the 
execution and three-view drawings. In the illustrations of the screen display, 
several occurrences of question marks preceded by numbers will be found. At 
these points in the program, data is required from the graphics tablet. These 
numbers correspond to locations on the three-view drawings from which the 
required data are taken. The numbers are composed of two parts. The first 
refers to a figure number while the second refers to a point on that figure, 
for example 5-A refers to point A in figure 5; a point labeled 4b-C would 
refer to point C on figure 4b. 

Figure 1 shows the first two steps in the execution of the program 
VORTAB. The first command gets procedure file PVORTAB 

GET, PV0RTAB/UN=21 4737C/ST=CPF 
The second command executes the procedure file. 

CALL, P VORTAB (NAME= EXAMPLE) 

EXAMPLE is the name of the file on which the resulting data will be stored. 

Any desired name can be used in place of EXAMPLE. 

Procedure File PVORTAB 

Procedure file PVORTAB, listed in Appendix A, is used to get PLANE, which 
is the binary form of program VORTAB, and also ensures that the files onto 
which VORTAB writes will be empty before execution. The binary program PLANE 
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is then executed. After the execution is completed, the three output files 
are saved as one packed file, under the name specified in the call statement. 
Then the two binary programs, READB and PLOTB (ref. 2), are obtained from 
permanent file and executed. This causes the configuration to be plotted 
allowing the adequacy of the modeling to be determined. 

Preparation and Use of the Three-View Drawing 

Use of program VORTAB (Appendix B) requires access to an interactive 
computer graphics terminal equipped with a graphics tablet. Before the pro- 
gram can be used, a number of items must be prepared. This program is designed 
for use with a configuration three-view, drawn to scale. The three required 
views are a top view, with the nose to the left, a side view, with the nose 
to the left, and a front view (see fig. 2). The drawing may be placed on the 
tablet in one piece or it may be separated into two or three separate sections. 
To obtain the greatest accuracy, the drawing should be as large as possible 
while fitting on the usable portion of the tablet. The precise orientation 
of the three views, or of each of the separate views, is unimportant; however, 
a line that would be horizontal on any of the views should Intersect the bottom 
of the tablet at an angle less than 45 degrees. The required configuration 
panels can be defined on either side of the aircraft X-Z plane, as symmetry 
is assumed, and all the panels will be stored in the VORLAX Input data files 
assuming the right-hand side of the configuration. Lines which are hidden in 
the front view will present problems in locating some panels. If the front 
and side views are aligned properly, the vertical locations of the hidden 
front-view panels can be taken from the side view. 

Organization of the VORTAB Program 

The VORTAB program is organized into six sections. The first section is 
used to determine the orientation of each view of the configuration three view. 
This section also determines the scaling factor to be used to convert the data 
from tablet units into aircraft units. The second through fifth sections are 
used to define individual major panels. These panels are divided Into four 
categories with one section of the program allocated for each. The categories 
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are; (1) wings and other flat panels that can have camber, (2) fuselage pan- 
els, (3) vertical flat panels, and (4) nacelle panels. As many panels as are 
needed to define the various configuration components may be defined with pro- 
gram VORTAB providing the limitations imposed by the VORLAX program are observed. 
The final section of the program is used to enter numerical data. These data 
have two purposes: to enter instructions used in processing the configuration, 

and to enter certain parameters that are used to define the configuration. 

After all of the data has been entered it is stored on disc in a permanent 
file. The name of the file is specified in the call to procedure file PVORTAB. 
After the data has been saved, two programs from PLOTIT (ref. 2) are executed 
which results in a plot of the aircraft configuration. This permits visual 
inspection of the resulting model. 

Initialization of the Program 

After the initial data and drawing preparation and before the program 
can be used to define the various configuration panels, a set of eight points 
must be entered to initialize the program. These points are used by the program 
to determine the location and orientation of each of the three configuration 
views. This section of the program requests the position of the two points 
called center leading point and center trailing point. These two points will 
always lie in the configuration plane of symmetry. If an odd number of fuse- 
lages are present, the center leading and trailing points are the nose and tail 
points of the central fuselage. If an even number of fuselages are present, 
the center leading point is the midpoint of a line joining the nose points of 
the two most central fuselages, and the center trailing point is the midpoint 
of a line joining the tail points of these same two fuselages. In a config- 
uration that has no fuselage, e.g., a flying wing, the center leading and 
trailing points are taken as the leading and trailing point of the wing center- 
line. 

The screen display which is generated during the initialization is shown 
in figure 3(a), and the curser locations are shown in figure 3(b). 
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First the orientation of the top-view is determined, by sending the 
locations of the central leading point (3b-A) and the central trailing point 
(3b— B) to the computer. The distance between these points in aircraft units 
Is then typed In at the keyboard. The front-view is Initialized next by 
specifying two points on a line parallel to the Y-axis, 3b-C and 3b-D, and the 
central leading point 3b-E, The sideview is initialized last, by specifying 
two points on a line parallel to the X-axis, 3b-F and 3b-G, and the central 
leading point 3b-H. 


Defining an Aircraft Configuration 

For this discussion, aircraft configurations will be divided into three 
catagories, aircraft with one fuselage, with an even number of fuselages, 
and with an odd number of fuselages greater than one. The following describes 
the order in which the panels are to be defined for each of the three cases. 

Single fuselage. - For a configuration with one fuselage, the panels are 
defined in the following order: wing panels (as many as necessary), other 

winglike or flat panels, fuselage panels, vertical panels and, finally, nacelles. 

Even number of fuselages. - If a configuration has an even number of 
fuselages, the first panels to be defined are the innermost wing panels followed 
by the innermost fuselage and continuing outward. After all of the wings and 
fuselages have been defined, the other winglike panels are defined, followed 
by vertical panels and nacelles. 

Odd number of fuselages. - When a configuration has an odd number of fuse- 
lages greater than one, the first panels to be defined are the wing panels 
between the central fuselage and the next outward fuselage. The central fuse- 
lage will then be defined, followed by the next outward fuselage. The other 
wing panels and fuselages are then defined continuing outward. After all of 
the wings and fuselages are defined, other winglike and flat panels can be 
defined followed by vertical panels and nacelles. 
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Description of Lifting Surfaces 

The first type of panel to be considered is the wing. Other panels such 
as canards and filler panels are special cases of wings. In the program these 
are called flat panels. If no panels of this type are needed, this section 
of the program can be bypassed. There is a separate section within the program 
to handle panels parallel to the plane of symmetry. 

The first question which appears in figure 4(a) asks if flat panels are 
to be defined. If so, the first panel is assumed to be the innermost wing 
panel . 

Location. - To define a flat panel, a set of six points is required from 
the graphics tablet as shown in figure 4(b). They are the root leading point 
(4b-A), root trailing point (4b-B), the vertical location of the topmost point 
of the panel root (4b-C), the panel tip leading point (4b-D), tip trailing 
point (4b-E), and the topmost point of the panel tip (4b-F). The required 
numerical data, entered through the keyboard, is shown in figure 4(a). 

Camber. - Horizontal panels can be either cambered or uncambered. If 
they are uncambered, this section of the program will be skipped. In the case 
of the wings, if the innermost wing panel is cambered all of the connected 
wing panels will be cambered. Distances from the leading edge for camber 
data can be defined in one of two ways, either as percent chord or linear 
measurement. If the data is entered in linear units, the program converts 
them to percent chord. In both cases the data is entered from the keyboard, 
and data for both the leading and trailing edges must be included. The 
camber values themselves are always input in units of percent chord. 

When the locations are defined in terms of percent chord as in figure 
4(c), the data must be entered in sets of three numbers, the percent chord 
station, the camber at the root-chord line, and the camber at the tip-chord 
line. Data must be included for stations at zero percent chord and at 100 
percent chord. Up to 50 stations may be used. After all of the camber data 
has been entered a set of any three numbers must be entered, the first of 
which must be negative, to terminate the definition of camber data. 
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When the locations are defined In terms of linear measurement, figure 
4(d), all of the data for the panel root Is entered first. The numbers are 
entered in pairs, the first of which is the distance from the leading edge 
while the second is the camber. The first pair must be for the leading 
edge while the last must be for the trailing edge. 

After the required data has been entered, two additional numbers must 
be transmitted, the first of which is negative. This instructs the program 
that no more data is to be entered for the panel root. Data is then entered 
for the panel tip, in the same manner as for the root. If a flat panel is 
to be uncambered, the display will appear as it does in figure 4(e). 

Error recovery. - Because the possibility exists that errors will be 
made while entering panel data, the last question displayed on the terminal 
display asks if the panel just defined is wanted. 

Additional flat panels. - After the first horizontal panel has been 
accepted, the Tektronix screen will appear as is shown in figure 4(f). Three 
options for the continuation of the program are offered at this point. The 
first, selected by entering an E, allows an extending panel to be added to 
the previously defined panel. If the previous panel was cambered, the 
extending panel must be cambered; if the previous panel was uncambered, the 
extending panel must be uncambered. The second option, selected by entering 
a N, makes it possible to define an entirely new panel. The final choice, 
selected by entering a Q, causes the program to preceed to the next type of 
panel. After each flat panel other than the first has been accepted or 
rejected, the display will appear as it does in figure 4(g). 

Description of Fuselage Panels 

Any number of fuselage panels can be defined within the limits imposed 
by the VORLAX program. Each of the fuselages is assumed to have an octagonal 
cross-section and may be composed of up to four major panels. 

The cross-section of the fuselage is initially assumed to be circular 
and is approximated by a regular octagon. This octagon is then modified so 
that the vertex nearest the wing root will be in contact with the wing root. 
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This allows fuselages with a noncircular cross-section to be approximated 
by an irregular octagon. The numerical data required for a fuselage is 
entered through the keyboard and is shown in figure 5(b). 

The first question, asked in figure 5(a), makes it possible te define 
an aircraft either with or without a fuselage. The first points entered In 
the definition of a fuselage are the nose point from the top-view, 5c-A, and 
from the side-view, 5c-B, as shown in figure 5(c). 

The next points, from the side-view, are entered in pairs, one from the 
top and one from the bottom of the fuselage with the bottom point directly 
below the top point. The pairs of points are shown as 5c-C-N and 5c-D-N, 
where N is the pair number. As many pairs as required can be entered within 
the limits of the VORLAX program. After the last desired pair of points has 
been entered, the program will again ask for a top point. At this point a Q 
should be entered which indicates that all of the point pairs have been entered. 
At this time it is possible to reenter the data If an error has been made. 

The desired number of fuselage panels is then specified up to a maximum of 
four. Each panel consists of an octagonal segment of a fuselage. The length 
of each segment, other than the one which ends at the tail point of the fuse- 
lage, is specified by entering its end points from the tablet. An example 
of such a point is 5c-E. 

The next data to be entered are the number of streamwise vortices. If 
an error has been made, the next question allows the panel to be rejected 
and the length and number of vortices to be redefined. If the fuselage was 
located at the center of the configuration the program will ask, IS ANOTHER 
FUSELAGE DESIRED? This would be the second fuselage of a multi-bodied 
configuration with an odd number of fuselages. 

Description of Vertical Panels 

A vertical panel is one which is parallel to the plane of symmetry, 
e.o., a fin or rudder. A vertical panel could also be defined in the section 
of the program used to define flat panels. The purpose of this section is to 
simplify the definition of vertical surfaces and avoid problems caused by 
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hidden lines. 

The questions which appear on the screen during definition of a 

vertical panel are shown in figure 6(a). The first question allows the user 

the option of either defining or not defining vertical panels. 

Four points are required from the side-view: the bottom leading point 

(6b-A), the bottom trailing point (6b-B), the top leading point (6b-C), and 
the top trailing point (6b-D). One additional point is required from the 
top view (6b— E ) and is used for spanwise location only. The numerical data 
which is required are the numbers of spanwise and streamwise vortices. At . 
this point the panel can be rejected if errors have been made while defining 

it. The final question, DO YOU WANT ANOTHER VERTICAL PANEL?, allows additional 

vertical panels to be defined. 

The fifth section of the program is used to define nacelles which are 
represented as uncambered hollow octagonal cylinders. Figure 7(a) shows the 
terminal display which results from the definition of a nacelle. The first 
question allows this portion of the program to be skipped. 

The first two points are taken from the top view and are the two front 
corners of the nacelle, points 7b-A and 7b-B. The next points, also from 
the top-view, are taken from one side of the nacelle and are labeled Cl 
through Cn. These points are used to determine the relative diameters of 
the nacelle in the streamwise direction. The last of these points must come 
from one of the two rear corners to determine the nacelle length. A Q is 
sent to the program to indicate that the last of these points have been entered. 
The final point, from the front-view, is the topmost point of the nacelle 
(7b-D). The number of streamwise vortices must be entered through the key- 
board. A nacelle can be rejected if an error is made in entering the data. 
Additional nacelles can be defined if necessary. 

Numerical Data 

This section of the program is used to enter parameters which define 
some aspects of the aircraft configuration and its operating conditions. It 
is arranged to correspond to cards in the VORLAX (ref. 1) input data deck 
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with the associated card numbers. Most of the parameters which must be 
entered In this section have names In parentheses which are the same as the 
names used In reference 1. The card number and variable name can be used 
to look up an explanation of the variables which Is more extensive than that 
which appears in this program. After all of the data for a card has been 
entered, the program asks if the data is to be reentered for that card. This 
allows erroneous data to be corrected. Data for any card can be reentered 
as many times as necessary. 

Figure 8(a) shows a possible screen display generated while entering 
data for the first six cards. Figures 8(b) and 8(c) show two possible 
screen displays generated while entering data for the last three VORLAX data 
cards. 

With the termination of this section of the program, execution of program 
VORTAB is complete. 

GRAPHIC OUTPUT 

At this point procedure file PVORTAB will execute programs READB and 
PLOTB which plot the aircraft configuration. 

Using figure 2 as a sample configuration, program VORTAB was used to 
generate the numerical data shown in figure 9. Utilizing the plotting rou- 
tines of reference 2, various views of the configuration have been plotted 
and are shown in figure 10. 
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/GET, PU0RTAB/ST-CPF/UN-214737C . 
/CALL (PUORTABC NAME ^EXAMPLE ) ) 


Figure 1. - Required two statements to retrieve and execute program VORTAB 
from mass storage. 
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Figure 2. - Relative positions of the required three views of a Supersonic Cruise Research (SCR) configuration 


CO 


FRON THE TOP VIEW 
- , „ CENTER LEADING POINT 

3 b-A ? 

_ , ^ CENTER TRAILING POINT 

3 b - B t 

INPUT THE LENGTH OF THE FUSELAGE (A HUNDER)* 

? 1M 

THREE POINTS FROfl THE FRONT UICU PRC ENTERED NEXT* 
FIRST POINT AND SECOND POINT REFER TO 8 POINTS ON 
A LINE PARALLEL TO THE V-AXIS. 


3b-C 
3 b-D 
3b- E 


FIRST POINT 
T 

SECOND POINT 
7 

CENTER LEADING POINT 

7 

THREE POINTS FRON THE SIDE UIEU ARC ENTERED NEXT 
FIRST POINT AND SECOND POINT REFER TO TWO POINTS ON 
A LINE PARALLEL TO THE X-AXIS* 


3b-F 

3b-G 

3b-H 


FIRST POINT 
? 

SECOND POINT 
T 

CENTER LEADING POINT 
T 


(a) Display of messages necessary in the initialization process. 
Figure 3. - Program VORTAB initialization process. 
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ME FLAT PANELS DESIRE* 

enter v or r for yes, to nuffsho 
’the first panel is a uinq attaches to the FUSLAOE. 


4b-A 
4b-B 
4b-C 
4 b-D 
4b-E 
4b-F 


ROOT LEABINO POINT (FRON THE TOP UIEU) 

ROOT TRAILING POINT (FRON THE TOP UIEU) 

ENTER UIHG/CANAAD EOOE POINT (FRON THE FRONT UIEU) 

OUTER LEASING POINT (FRON THE TOP UIEU) 

TIP TRAILING POINT (FRON THE TOP UIEU) 

ENTER UING/CAHARS ESGE POINT (FRON THE FRONT UIEU) 

? NUN*ER OF SPAHUISE MORTICES (NUOR) 

? NUNBER OF STREAHU1SE MORTICES (RHCM) 

^LEADING ESGE SUCTION HULT1PLIER tSPC) 0.0<«SPC<»1.0, 
SPCCO.O FOR POLHANUS ANALOGY 

? ENTER THE ANGLE OF THE ROOT CHOPS LINE UN DEG. ) 

A POSITIVE ANGLE CAUSES A DEFLECTION DOUN. 

T EHTER THE ANGLE OF THE TIP CHORD LI« UN DEG.) 

T 0 


OR 


(a) Display of messages necessary in defining a flat panel 
Figure 4. - Definition of a wing. 
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I 

i 


DO YOU WANT CARDER IN THE PANELS 
ENTER Y OR T FOR YES, OR N OR F FOR NO 
7 V 

CAUSER IS ENTERED IN R-CHORD OR LINEAR MEASURE. 

DO YOU UANT CAUSER IN R-CHORD 

ENTER Y OR T FOR YES, OR N OR F FOR NO 

7 Y 

INPUT CANKER SATA BELOW. 

INPUT PERCENT CHORD < #.# TO STOP 
R CHORD. ROOT CAMBER. A TIP CAMBER 

? S*CHORD. ROOT CAMBER, A TIP CAMBER 

7 S. 1.4.3. 3.1 

R CHORD. ROOT CAMBER. A TIP CAMBER 
7 lt.«.6.*,S.« 

R CHORD. ROOT CAMBER, A TIP CAMBER 

7 2A.S.4 

R CHORD, ROOT CAMBER. A TIP CAMBER 

7 44 , 4 , 3 . 

R CHORD. ROOT CAMBER. A TIP CAMBER 
7 7R.2.1 

R CHORD. ROOT CAMBER. A TIP CAMBER 
7 1M, 1,1 

R CHORD. ROOT CAMBER, A TIP CAMBER 
7 -1,2.3 

DO YOU UANT THE PREVIOUS PANEL 
ENTER Y OR T FOR YES, OR N OR F FOR NO 
7 V 


(c) Display of messages necessary in defining panel camber in percent chord. 


Figure 4 


Continued 



DO YOU UAHT CANBCff IN THE PANELS 

ENTER V OR T FOR YES, OR N OR F FOR NO 

? V 

CANBER IS ENTERED IN X-CHORO OR LINEAR NEASURC. 

DO YOU WANT CARDER IN S-CHORD 

ENTER V OR T FOR YES, OR N OR F FOR NO 

? N 

INPUT CARDER DATA DELOU 

TO STOP INPUT DISTANCE FROA LEADING EDGE < M 
ENTER DIST. FROR LEAD. EDGE, AND ROOT CARDER. 

7 0 0 

ENTER DIST. FROfl LEAD# EDGE# AND ROOT C MISER. 

? 4.3,3 

ENTER DIST. FROfl LEAD. EDGE# AND ROOT CARDER. 

? 1B#3 

ENTER DIST. FROfl LEAD. EDGE# AND ROOT CAflBER. 

? 30 , 1.6 

ENTER DIST. FROfl LEAD. EDGE# AND ROOT CAflBER. 


? 50#. 1 

ENTER DIST. FROfl LEAD. 


EDGE# AND ROOT CAfltCR. 


? -I#l 

STATION NURBER 1 

ENTER DIST. FROfl LEAD. EDGE 4 TIP CAflBER. 

? •#• 

STATION NURBER 3 

ENTER BIST. FROfl LEAD. EDGE 4 TIP CAflBER* 
? 2.15*2.3 

STATION NURBER 3 

ENTER DIST. FROfl LEAD# EDGE 4 TIP CARDER. 
? 5,2 

STATION NURBER 4 

ENTER DIST. FROfl LEAD. EDGE 4 TIP CARDER. 
? IS# .75 

STATION NURBER S 

ENTER DIST. FROfl LEAD. EDGE 4 TIP CAflBER. 
? 2S#.B 

DO VOU WANT THE PREVIOUS PANEL 
ENTER V OR T FOR YES, OR N OR F FOR NO 
? N 


(d) Display of messages necessary in defining panel camber in linear measurements. 

Figure 4. - Continued. 
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* «VOU UftNT C MOCK IN THE PANELS 

enter v on t rot yes, on n or f for no 

PO YOU UAHT TIC PREVIOUS PANEL 
ENTER Y OR T POR YES* OR N Oft F FOR NO 
IT 


(e) Display of messages necessary if a flat panel is not cambered. 


ONE Or THREE POSS1SILIT1ES RUST K CHOSEN NEXT. 

(1) ft PANEL CAN K ADDED TO THE OUTER EPOC OF THE 
PREVIOUS ONE C I.E. ft NVLT I -PftNELCO UINC), 

(2) ftN ENTIRELY NEU FiftT PftNEL CRN >E ENT ERED, 

<3 > THE NEXT TYPE OF PftNEL CftH SC ENTERED. 

FOR ftN ENTCNDINQ FLftT PftNEL INPUT. . •£ (EXTENDXHQ). 

FOR ft NEU FLftT PftNEL INPUT N (NEU PftNEL >. 

TO STOP FLftT PftNCLS INPUT ..0 (QUIT). 


(f) Display of choices after the first horizontal panel has been accepted. 


FOR ftN CNTENSINQ FLftT PANEL INPUT 

FOR ft NEU FLftT PftNEL ItTUT 

TO STOP FLftT PANELS INPUT 

? O 


E 

N 

.0 


(EXTENDI NO). 
(NEU PANEL). 
(SUIT). 


(g) Display of choices for subsequent horizontal panels. 

Figure 4. - Concluded. 
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DO YOU UftNT * FUSELAGE 

ENTER V OR T FOR YES, OR N OR F FOR NO 

7 Y 


ENTER NOSE POINT FRON TOP VIEW 

5c-A i 


(a) Display of message necessary in defining a fuselage 


Figure 5. - Definition of a fuselage. 



r>o 

ro 


TO NMCC 0 FUSCLMC SOSV CNTCN MINT* FNON TMf SIM UIIU . TMJE I NOSC JOINT IS nilST. TMt MWjT MINTS MI CNTCRO IN 
PNINS, TNI FIOTT FNON THK TOO Of TNI FUSCLMI* I NI **COM J[SM TM WTTOW SIMCTtV SCLOW TNI FINNT. 

TW LMT TWO COM FNON TNI IH# OF TNI FUSCLMt. FFTW TWSt TWO CNTCN 0 TO DIO 


5c-B 


NOK F01NT 


|C- 

>c- 


V “ • I 

8=il 


TOF POINT 

5C-C-1 f SOTTON POINT 
" TOP POINT 
S OTTON POINT 
TOP POINT 
t SOTTON POINT 
,T TOP POINT 
T SOTTON POINT 
T TOP POINT 
f SOTTON POINT 
f TOP POINT 
T S OTTON POINT 
T TOP POINT 
f SOTTON POINT 
T TOP POINT 
T SOTTON POINT 
T TOP POINT 
T SOTTON POINT 
T TOP POINT 
T 
T 


5c-C-9 

5C-D-9 


SO VOU WONT TO NC-CNTCN THIS ONTO 

OfTCR VMTFM VIS* ON N ON F FON NO 

VuNSCN OF FUSCLA6C PMCLS 

T INTtt SUM POINT FON PNNtt S 1 , (FNON TNI TOP OICU) 


5c-E t 


NUNKN OF CNONSWISt UONTICICS (NNCV) 

FON PMCL NUHSCN t 
T II 

SO VOU WNNT TNI PMWIOUS PMCL 
CNTIR V ON T FON VIS, ON N ON F FON NO 
t V 

NUHSCN OF CNONMXK UONTICICS (NNCU) 

FON PMCL NUHSCN I 
T II 

00 VOU WNNT TME PNCVIOUS PMCL 
orrat V ON T FON VIS* ON N ON F FON NO 
T V 

IS MWTNCN FUSCLAOI NCOUINCS 
CNTCN V ON T FON VIS. ON N ON F FON NO 
T V 


(b) Display of messages necessary in defining a two panel fuselage. 


Figure 5. - Continued. 
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PQ you uant vertical panel* 

opt* V M T FOR «l, OR N OR F FOR NO 

T Y 

g b _ A ROTTON LEADIN8 POINT (FRON THE SIK VIEW) 

6 b - B ?*”* TRAILIHO POINT (FRON THE SIK VIEW) 

g b _ c TOP LEADINO POINT (FRON THE SIK VIEW) 

6 b - D T* T ** ILI,W MIHT <«WH THE SIK VIEW) 

6b-E? P0IKT F0R LOCATION I FRON THE TOP VIEW) 

^NURKR OF SPNHVISE VORTICES (NVORI 

^CMNKR OF STREAHUISE VORTICES <RNCV> 

PO VOU UANT THE PREVIOUS PNHEl 
^EHTER V OR T FOR YES, OR N OR F FOR NO 

RO VOU UANT ANOTHER VERTICAL PANEL 
^ENTER V OR T FOR YES. OR N OR F FOR NO 


(a) Display of messages necessary in defining vertical panels 
Figure 6. - Definition of a vertical surface. 




7b-A 

7b-B 


7 b-C 
7b-D 


M VOU WANT NACELLES 

ENTER V OR T rot VES, Ot N OR F FOR NO 

t Y 

mo* THE TO* VIEW 

INNER OR OUTER FRONT NACELLE POINT 

? 

RENAXHINQ FRONT NACELLE POINT 
? 

ENTER SODV POINTS ALONQ ONE SIDE (EXCLUDING FRONT 
POINTS!. PROQRCSSXNQ FROM FRONT TO REAR AFTER THE 
LAST POINT HAS SEEN PROCESSES ENTER 0 FOLLOWER SV 
A CARR I ACE RETURN* 

M T 0 

TOP NACELLE POINT CFRON THE FRONT VIEW) 

? 

NUNBER OF CHORSWISE VORTICES CRNCU) 

T 9 

SO VOU WANT THE PREVIOUS PANEL 
ENTER V OR T FOR YES, OR N OR F FOR NO 
? V 

IS ANOTHER NACELLE SESIRES 
ENTER V OR T FOR YES, OR N OR F FOR NO 
V N 


(a) Display of messages necessary in defining a nacelle 

. - Definition of a nacelle. 


Figure 7 




no 

co 


THE CARD NUMBERS CORRESPOND TO THE CARO NUMBERS 
AS GIVEN IN THE DOCUNENT POM THE VOBLAX PROGRAM. 
ENTER THE DATA FOR CARD NUMBER 1. 

ENTER TITLE FOR THE AIRCRAFT CONFIGURATION 
? THIS IS A SAMPLE EXECUTION 
DO VOU UANT TO RE-ENTER THE DATA FOR CARO 1 
ENTER V OR T FOR YES, OR N OR F FOR HO, 

? T 

ENTER THE DATA FOR CARO NUMBER I. 

ENTER TITLE FOR THE AIRCRAFT CONFIGURATION 
? THIS IS STILL A SAMPLE EXECUTION. 

DO VOU UANT TO RE-ENTER THE DATA FOR CARO 1 
ENTER V OR T FOR YES, OR N OR F FOR NO, 

? N 

ENTER THE DATA FOR CARD NUMBER 3. 

METHOD OF SOLVING THE SYSTEM OF EQUATIONS RELATINO 
BOUMDRV CONDITIONS TO UORTICITY STRENGTH ( I SOLVE ) 

• FOR GAUSS-SEXDEL 

I FOR PURCELL'S VECTOR ORTHOGONAL! ZAT I OH 

? CHORBUISC SPACING OF VORTICES (LAX) 

• FOR COSINE LAU SPACING 
1 FOR LINEAR SPACING 

7 1 

SPANUISC SPACING OF VORTICES (LAV) 

• FOR COSINE LAU SPACING 
I FOR LINEAR SPACING 

7 0 

INPUT THE OVER-RELAXATION PARAMETER (REXPAR) 

INPUT M FOR MACHINE DETERMINED OPTIMUM 
*»• <- REXPAR < ,00 
7,1# 

HEIGHT ABOVE GROUND (HAG), 

ENTER M FOR NO GROUND EFFECT 

7 • 

LONGITUDINAL VORTEX UAICE FLOTATION FACTOR (FLOAT* ) 
7 ,75 

LATERAL VORTEX UACE FLOTATION FACTOR (FLOATY) 

T .GJ 

MAXIMUM NUMBER OF ITERATIONS FOR THE GAUSS-SEXDEL 
RELAXATION METHOD (XTRNAX) 

IF t IS ENTERED DEFAULT OF 99 UILL BE USED 
7 3 

DO YOU UANT TO RE-ENTER THE DATA FOR CARO 2 
ENTER V OR T FOR YES, OR N OR F FOR NO, 

7 H 

ENTER THE DATA FOR CARD NUMBER 3. 

NUMBER OF MACH NUMBERS TO BE ANALYZED (NMACH) 

<-7 


INPUT MACH NUMBER t 1 

7 »o'vOU UANT TO RE-ENTER TNE BATA FOR CANS 3 
CNTCR Y OR T FOR VW, ON H ON F FON NO. 

7 N 

CNTCR THC DATA FON CARD KUMCN 4. 

NUNBCN OF ANOLES OF ATTACK (NALPHA) 

<•7 

7 JNFUT ANGLE * 1 IN SCONCES 

? IHFUT ANGLE ON IN OCONCCS 

7 S0 VOU UANT TO RE-ENTER THC DATA FON CANS 4 

ENTER V ON T FON YES, ON N ON F FON NO. 

T N 

ENTER THE DATA FOR CARD NUHDCR 5. 

A9VNNCTRIC FLIGHT ON CONFIGURATION FLAG (LATNAL) 

• FON SYMMETRIC FLIGHT AND CONFICUNATION 
I FON ASYMMETRIC FLIGHT AND/ON COKFIQUNATION 
7 t 

SIDESLIP ANGLE IN DEGREES (PSD 
7 8 

PITCH RATE IN DEQRCES/SECOND (PITCHO) 

? NOLL RATE IN DCOREES/SCCONB (KOLLO) 

7 4 

VAU RATE IN DCGNCCS/SCCOND (VAUG> 

7 S 

REFERENCE FREE STREAM VELOCITY (UINF) 

IF S.t IS ENTERED DEFAULT OF t.S IS USED 

7 D0 VOU UANT TO RE-ENTER THC DATA FON CARD 6 
ENTER V ON T FON VIS. OR N OR F FON NO. 

7 H 

ENTER THE DATA FOR CARD KUNSER 6. 

REFERENCE MEA FON FORCE AND FOMENT COEFFICIENTS (SNEF) 

7 piTCHIN0 MOMENT COEFFICIENT REFERENCE LENGTH 
USUALLY MEAN AERODYNAMIC CHORD LENGTH (CM*) 

7 S 

X-COCNDINATC OF MOMENT REFERENCE POINT (XDAN) 

^-COORDINATE OF MOMENT REFERENCE POINT (ZDAR> 

? TOTAL UINC SPAN IN UNITS OF SNEF AND CDAR (USPAH) 

7 D0 # V0U UANT TO RE-ENTER THE DATA FOR CARD • 

ENTER Y OR T FOR YES. OR N OR F FOR NO. 

7 N 


(a) Display of messages necessary in defining the first set of numerical data. 
Figure 8. - Completion of configuration definition and flight conditions. 



THE DATA FOR THE LAST 3 CMOS XS INPUT NEXT, 

THESE CMOS ME NUMBERED 19, tt, MO 21. 

ENTER DATA FOR CARO MUR9ER 19. 

NUN9ER OF X- STATIONS TO DEFINE SPATIAL FLOW FIELD 
SURVEY GRID < 21, (NXS) 

? 2 

NUMBER OF Y-STATXCHS TO DEFINE THE SUIT LINE 
VALUES OF THE SURUEV GRID <NVS> 

V 2 

NUMBER OF 2-STATIONS TO DEFINE THE HATER LINE 
VALUES OF THE SURUEV GRID (N2S) 

7 2 

DO YOU UANT TO RE-ENTER THE DATA FOR CARD 19 
ENTER V OR T FOR VCS, OR N OR F FOR NO, 

? N 

ENTER THE DATA FOR CARD HUMBER 99 
X -STATION VALUES FOR SPATIAL FLOW FIELD (XS) 
i ENTER X -STATION NUNKR 1 
? 1 

ENTER X-5TATIOM NUNSER 2 

? 2.2 

DO YOU WANT TO RE-ENTER THE DATA FOR CMD 29 
ENTER V OR T FOR VES, OR N OR F FOR NO, 

? N 

ENTER DATA FOR CMD NUNKR 21 

BEGINNING OF GRIB IN BUTT LINE DIRECTION (VNOT) 

T 3 

V-SPACING OF THE GRID (PELT AY) 

? 2 

KG INNING OF QRID IN WATER LINE DIRECTION (ZNOT) 

! T 3 

S 2-SPACING OF THE GRID (DELTAZ) 

7 ^ 

10 YOU WANT TO RE-ENTER THE DATA FOR CARD 21 
ENTER Y OR T FOR YES, OR N OR F FOR NO, 

7 N 

EXECUTION OF THE PROGRAM IS COMPLETED 


(b) Display of messages necessary in defining the last three cards of program 
VORTAB when the first card is non-zero. 


Tt* MTU FO* THE LAST 3 CMN IS INPUT NEXT, 
THESE CMOS ME NUH K R E S IS, W. MB II. 


enter mta for cmo number is. 

NUNKR OF X-STATIOMS TO DEFINE SPATIAL FLOW FIELD 
SURUEV QRID < 11, <HXS> 

V 0 

SO VOU WANT TO ME -ENTER TIC SATA FOR CARO 10 
ENTER V OR T FOR VES, OR N ON F FOR NO, 

T N 

EXECUTION OF THE PROQNA H IS CONFUTES 


(c) Display of messages necessary in defining the last three cards of program 
VORTAB when the first card is zero. 

Figure 8. - Concluded. 
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Figure 9. - Numerical definition of a Supersonic Cruise Research (SCR) configuration 
presented in figure 1. 
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Figure 9. - Continued. 
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Figure 9. - Concluded. 



(a) Front view. 


Figure 10. - A VORTAB representation of a SCR configuration. 



(b) Side view. 


Figure 10. - Continued. 






APPENDIX A 

PROCEDURE FILE PVORTAB 


GET, PLANE/ST-CPF/UN-214737C . 
RETURN , TAPE3 , TAPE4 , NAME . 

PLANE . 

REPLACE, TAPE7-NANE. 

APPEND , NAME , TAPE3 , TAPE4 ) 

GET, NAME. 

PACK, NAME. 

REPLACE, NAME, 
rr fap 

GET, TAPEl -NAME. 

GET , READB , PL0TB/ST«CPF/UN»214737C 
READB . 

PLOTB. 

EXIT. 



APPENDIX B 


SOURCE LISTING OF PROGRAM VORTAB 


PROGRAM VORTAB ( I NPUT > OUT PUT # T APe 1 » T N PUT » T A P E2 # T APE 3> T AP E 4 > T AP E7 > 
C *♦ MAIN PROGRAM# CALLS SUBROUTINES THAT GENERATE PANELS. 

REAL MLINE# MS IDE 

COMMON /FUSS/XO, YO » FS CORD » F S L QP E » B FUS > NP AN S > FS C A LE 

COMMON /FRG/MLINE# BLINE# XXO# YYO 

COMMON /S IDE/MS I DE » BS I DE # XX# ZZ 

COMMON /ZXPR/ ZC (2#150) #XAF ( 150) >PHT (9) »R0 (9) 

DATA NPANS# YL0C/0#0.C/ 

C ** INITIALIZING ROUTINES ARE CALLED. 

CALL 1NITT(120) 

CALL TABINT ( 1 # 0# 0 ) 

CALL TERM ( 3# 4096 ) 

CALL CHRSIZC4) 

CALL FUSPRMS 
CALL ANMODE 

PRINT *# "THR E E POINTS FROM THE FRONT VIEW ARE ENTERED NEXT." 
PRINT *# "FIRST POINT AND SECOND POINT REFER TO 2 POINTS ON" 

PRINT *# " A LINE PARALLEL TO THE Y-AXIS." 

PRINT *#" " 

PRINT *, "FIRST POINT" 

CALL FNDPLAC(MLINE#BLINE#XXO> YYO) 

• CALL NEWLIN2 
CALL ANMODE 

PRINT *# "THREE POINTS FROM THE STOF VIEW ARE ENTERED NEXT" 

PRINT +# "FIRST POINT AND SECOND POINT REFER TO TWO POINTS ON" 
PRINT *#"A LINE PARALLEL TO THF X-AXIS." 

PRINT *t « " 

PRINT *# "FIRST POINT" 

CALL FNDPLAC(MSIDE#BSIDE#XX,ZZ) 

C ** ROUTINES FOR ENTERING VARIOUS PANELS ARE CALLED. 

1 CALL HTAIL (YLOC ) 

CALL NEWPAG 
CALL ANMODE 


VORTAB 

1 

VORTAB 

2 

VORTAB 

3 

VORTAB 

4 

VORTAB 

5 

VORTAB 

6 

VORTAB 

7 

VORTAB 

8 

VORTAB 

9 

VORTAB 

10 

VORTAB 

11 

VORTAB 

12 

VDRTAB 

13 

VORTAB 

14 

VORTAB 

15 

VORTAB 

16 

VORTAB 

17 

VORTAB 

18 

VORTAB 

19 

VORTAB 

20 

VORTAB 

21 

VORTAB 

22 

VORTAB 

23 

VORTAB 

24 

VORTAB 

25 

VORTAB 

26 

VORTAB 

27 

VORTAB 

28 

VORTAB 

29 

VORTAB 

30 

VORTAB 

31 

VORTAB 

32 

VORTAB 

33 



Co 

vo 



APPENDIX B. - Continued 


PRINT 

♦ >"DO YOU WANT A FUSELAGE'’ 

VORTAB 

3A 

CALL 

YESNOF(ISLT) 

VORTAB 

35 

PRINT 

* , « II 

VORTAB 

36 

IHISLT.LT. 0)G0 TO 2 

VORTAB 

37 

CALL 

FUSLAGEt YLOC ) 

VORTAB 

38 

lF(YLCC.EQ.O.O)GO TO 2 

VORTAB 

39 

GO TO 

1 

VORTAB 

AO 

2 CALL 

VTAIL 

VORTAB 

A 1 

CALL 

NACELLE 

VORTAB 

A2 

CALL 

NEWPAG 

VORTAB 

A3 

CALL 

ANMODE 

VORTAB 

AA 

CALL 

FIRS T6 (NPANS ) 

VORTAB 

A5 

CALL 

NEWPAG 

VORTAB 

A6 

CALL 

ANMODE 

VORTAB 

A7 

CALL 

LAST 3 

VORTAB 

A8 

CALL 

S YS TEM ( 0* " ’•) 

VORTAB 

A9 

CALL 

NEWPAG 

VORTAB 

50 

CALL 

ANMODE 

VORTAB 

51 

CALL 

FINITT (0* 3000 ) 

VORTAB 

52 

STOP 


VORTAB 

53 

END 


VORTAB 

5 A 



VORTAB 

55 


SUBROUTINE NEWLIN2 

VORTAB 

56 

** ADVANCES PRINTER TWO LINES. 

VORTAB 

57 

CALL NEWLIN 

VORTAB 

58 

CALL LINEF 

VORTAB 

59 

RETURN 

VORTAB 

60 

END 

VORTAB 

61 


VORTAB 

62 

SUBROUTINE RONE PNT < R J, RK ) 

VORTAB 

63 

** GETS A PAIR OF REAL COORDINATES FROM THE GRAPHICS TABLET. 

VORTAB 

6A 



APPENDIX B. 


COMMON /TABLET/ KTBMAX,KTBSND#« TRSFF 
DIMENSION IOUT13) 

DATA I OUT ( 1 )» IOUT (2 )/27> 33/ 
IGUT(3)«32+KTBMAX+KTBSND*KTBSEE 
1 CALL TOUTS T ( 3> I OUT ) 

CALL GETPNT ( l» J>K) 

IF( I.GT.3DGO TO 1 

IOUT( 3 )-60 

CALL T0UTST(3, IOUT) 

RJ«FLOAT( J) 

RK-FLOATIK ) 

RETURN 

END 


SUBROUTINE FRNTLQC ( Z1>ARAY, IS 17) 

C ** GETS A POINT AND LOCATES IT XJTH RFSPFCT 
C ** FRONT VIEW. 

REAL MLINE 

COMMON /FRO/MLINE>BLINE»XXO, YYO 
DIMENSION ARAY(ISIZ) 

DIMENSION AO ( 3 ) 

DATA AO/1 OH (FROM THE > 10HFR0NT VIFW.3H) 

CALL NEWLIN2 

CALL AOUTST(ISIZ,ARAY) 

CALL A0UTST(23> AO) 

CALL RONE PN T ( X> Y ) 
IF(ABS(MLINE).LT..C01)G0 TO 5 
B * Y + X /ML I NE 

XL INE * ( B-BL INE )/ (MLINE + 1. O/ML IN E ) 

YLINE-MLINE+XLINE+BLINE 

Z1«SCRT((XLINE-X)**2+(YL1NE-Y)**2) 

IF(Y.LT.YLINE)Z1*-Z1 

RETURN 
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VORTAB 76 

VORTAB 77 

VORTAB 78 

VORTAB 79 

TO THE VORTAB 80 

VORTAB 81 

VORTAB 82 

. VORTAB 83 

VORTAB 84 

VORTAB 85 

/ VORTAB 86 

VORTAB 87 
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-p» 

ro 

5 Z1-Y-YY0 
RETURN 
END 
C 


SUBROUTINE LQCATE(X*Y) 

C ♦ ♦ LOCATES POINTS WITH RESPECT TO FUSELAGE CENTER LINE* 
C ♦♦FROM THE TOP VIEW . 

COMMON /FUSS/XO* YO* FSCORD* FS LOPE » BFUS *NP ANS * FSC ALE 
SIGN-1.0 

I F ( ABS (FSLOPE ) .LT. .001 )G0 TO 5 
B2-Y+X/FSL0PE 

XL I NE = ( B2 rB FUS ) / ( F SLOP E +1 .O/FSLOP E ) 

YLINE-FSLOPE*XLIN£+BFUS 

IF (Y.LT. YLINE ) SIGN— SIGN 

Y-SQRT ( ( XLINE-X )♦♦£+( YL I NE-Y ) +*2 ) *STGN 

X -S QRT ( ( YLINE-Y0>*^2+ ( XL INE-XO)**?) 

RETURN 
5 X-X-XO 
Y-Y-YO 
RETURN 
END 
C 
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SUBROUTINE LOCATES (X*Z ) 

C ♦♦ LOCATES POINTS WITH RESPECT TO THF SIDE VlFW. 
REAL MSIDE 

COMMON/SIDE/MSIDE* BSIDE*XX,ZZ 
IFtABS(MSIDE) .LT..001JG0 TO 2 
SIGN-1.0 
B2-Z+X/MSIDE 

XLINE*(B2-BSIDE)/(MSIDE+1.0/MSIDF) 

ZLINE-MS IDE+XLINE+BSIDE 
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IF(Z.LT.ZLINE)SIGN*-1 .0 

VORTAB 

129 

Z»SQRT((XLINE-X)**2+(ZLINE-Z)**?)*STGN 

VORTAB 

130 

X-SORT( ( ZLINE-ZZ)**2 + <XLINE-XX)**2 ) 

VORTAB 

131 

RETURN 

VORTAB 

132 

2 X«X-XX 

VORTAB 

133 

Z-Z-ZZ 

VORTAB 

1 34 

RETURN 

VORTAB 

135 

END 

VORTAB 
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SUBROUTINE R ES C ALE ( X> Y, Z > CD, X 1» Y1 , 7 1 , CD1 » P D L, N VOR ) 

VORTAB 

138 

*♦ RE-SCALES THE OATA FROM SCREEN UNITS TO USER UNITS. 

VORTAB 

139 

COMMON /fuss/xo»yo>fscord*fslupf»rfus»npans»fscalE 

VORTAB 

1 90 

COMMON /ZXPR/ ZC12>150)>XAF(150)»PHI(9)»R0(9) 

VORTAB 

191 

X«X*FSCALE 

VORTAB 

192 

Y«ABS (Y+FSCALE) 

VORTAB 

193 

Z«Z*FSCALE 

VORTAB 

199 

CD*CD*FSCALE 

VORTAB 

195 

X1»X1*FSCALE 

VORTAB 

196 

Yl-ABS (Yl + FSCALE ) 

VORTAB 

197 

Z1*Z1*FSCALE 

VORTAB 

1 98 

CD1«CD1*FSCALE 

VORTAB 

199 

IF (PDL.EQ .0.0 JRETURN 

VORTAB 

150 

DO 10 J ■ 1 > N VOR 

VORTAB 

1 51 

10 RO(J)=RG(J)*FSCALE 

VORTAB 

152 

RETURN 
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END 
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SUBROUTINE LINFEED(N) 

*+ ADVANCES PRINTER N LINES. 
IF (N.EQ.O)RETURN 
CALL CARTN 
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DO I J * 1 > N 
1 CALL LINE F 
RETURN 
END 
C 


SUBROUTINE FNOPLAC (ML INfc,BLINE»X,Y) 

C ** DETERMINES SLOPE OF REFERENCE LIN? AND LOCATES 
C ** NOSE POINT FOR THE FRONT AND SIDF VTFWS. 

REAL MLINE 

DIMENSION A 1 ( 2 ) , A2 ( 2 ) 

DATA Al/IOHSECOND P0I,10HNT / 

DATA A2 / 10HCENT ER LEA* 10HDING POINT/ 

CALL SYSTEM 
CALL LINFEED(A) 

CALL R0NEPNT(X1,Y1> 

CALL NEWLIN2 
CALL AOUTS T (20, A1 ) 

CALL RONEPNT ( X2, Y2 ) 

MLINE«(Y2-Y1)/(X2-X1) 

CALL NEWLIN2 
CALL AOUTS T (20, A2 ) 

CALL RONE PNT ( X, Y ) 

BLINE-Y-MLINE+X 

RETURN 

END 

C 


SUBROUTINE FUSPRMS 

C ** GETS FUSELAGE LENGTH AND SLOPE TO BE USED LATER* 

C ** AND NOSE POINT FROM THE TOP VIEW. 

COMMON /FUSS/XO,YO,FSCGRD,FSLOPE,RFUS,NPANS,FSCALE 
DIMENSION N0SEPNT<2), TAILPNT(3) ,TOPVTEW( ?) 


VORTAB 

160 

VORTA8 

161 

VORTAB 

162 

VORTAB 

163 

VORTAB 

16* 


VORTAB 

165 

VORTAB 

166 

VORTAB 

167 

VORTAB 

168 

VORTAB 

169 

VORTAB 

170 

VORTAB 

171 

VORTAB 

172 

VORTAB 

173 

VORTAB 

17* 

VORTAB 

175 

VORTAB 

176 

VORTAB 

177 
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DATA N0SEPNT/10HCENTER LEA* 10HDTNG PQTNT/ 

DATA TAILPNT/IOHCENTER TRA# IOHIL ING P0IN,1HT/ 

DATA TOPVIEW/IOHFRCM THE T,7HOP VIFW/ 

C ** NOSE POINT. 

CALL NEWPAG 

CALL AOUTST (17»T0PVI£W) 

CALL NEWLIN 

CALL A0UTST(20>N0SEPNT) 

CALL RONEPNT ( XO » YO ) 

C *♦ TAIL POINT. 

CALL NEWLIN2 

CALL AOUTST (21>TAILPNT) 

CALL RONEPNT(X,Y) 

FSLOPE«(Y-YO)/(X-XO) 

BFUS*Y-FSLOPE*X 

C *♦ CALCULATES THE LENGTH OF THE FUSELAGE IN RASTER UNITS. 
FSCORD-SQRT( ( X-XO ) **2 + ( Y-YO ) **?. ) 

I F ( X . LT. lOOOO. >G0 TO 5 
X-1.0 

5 PRINT " 

PRINT *> " " 

PRINT *> " IN PUT THE LENGTH OF THF FUSELAGE (A NUMBER).” 

10 READ ( 1> * ) RLENGTH 
IF ( EOF (II ) 10# 11 

11 IF (RLENGTH. GT .-10. )G0 TO 6 
STOP 7 

C ♦* CALCULATES SCALE FACTOR TO CONVERT FROM PASTER 
C ♦* TO USER UNITS. 

6 FSCALE-RLENGTH/FSCORD 
CALL NEWLIN2 

CALL NEWLIN2 

RETURN 

END 
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197 
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199 
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201 
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207 
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210 
211 
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SUBROUTINE HTAIL(YLOC) 

C ** THIS INPUTS NGN-VERTICAL FLAT PANELS. 

COMMON /FUSS / XO / YO / FSCCiRD/FSLOPP » RF'JS »N PAMS /FSC ALE 
COMMON /WING/ WX ( 3 ) /WY ( 3 ) zWZ ( 3) » MOP p 
COMMON /Z XPR / ZC (2/150)/XAF (150), PHI (0) ,R0 (9) 
DIMENSION IANS ( 15 ) / AO (2 )z AI (2)/ A? <2 ) . AIR (2 ) 

DIMENSION F TP (3 ) . TOP (2 ) 

REAL NVOR 
LOGICAL FIRST 

DATA PIE / ITS/IQUANT/ISYNT/NPP/PDL /3 .1415926/4*0/0.0/ 
DATA FIRST/ IANSERz ISLT2/N02/ . TRUE. »72»1/-1 / 

DATA AO / lOHROOT LEADI/10HNG POINT / 

DATA T0P/1CH ( FROM THE /10HT0P VTFW) / 

DATA A1/10HRG0T TRAIL/10HING POTNT / 

DATA A1B/10HTIP TRAIL 1/ 9HNG POINT / 

DATA A2/10HOUTER LEAD/ 10HING POINT / 

DATA FTP/10HENTER W I NG/ 10H /C ANARO F0/10HGE POINT / 
DATA IP AN EL/1/ 

C ** SOME PARAMETERS ARE INITIALIZED. 

85 I P AN 3 * 1 
I PAN4=-1 
I PAN2*1 
CALL NEWPAG 
CALL ANMODE 
3 CALL NEWLIN2 

C ** DETERMINES IF FLAT PANELS ARE DESIRED. 

IF ( I PANEL . EQ. 1 ) GO TO 1 
IF ( YLOC.EO.O.OIGO TO 1 

8 PRINT */ "DO YOU WANT THIS FLAT PANEL ATTACHED" 

CALL YESNOF (IPAN3) 

IF(N02.LT.0)CALL L INF EED ( 3 ) 

IPAN4-1 

I F ( I APN3 ) 9/ 9/ 11 

1 PRINT */ "ARE FLAT PANELS DESIRED" 

CALL YESNOF (MODE ) 
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. - Continued 


c 


IMMODE.L 
& CALL NEWL 
** DETERMINES 
PRINT *,'« 
CALL SYST 
10 IF ( IPANA ) 
9 ITESTW-0 


T . 0 ) RETURN 
IN2 

IF WING AND FUSELAGE ARP TO BF CONNECTED. 

THE FIRST PANEL IS A WING ATTACHED TO THE FUSLAGE. 

EM (0, " »> 

9,9,8 


it 


NAP-0 

IANSR2-IANSER 
11 CALL NEWLIN2 
C ** GETS ROOT LEADING POINT. 

CALL AOUTST ( 19, AO) 

CALL AOUTSTt 19, TOP ) 

CALL RONE PNT ( Xl , Y1 ) 

CALL LOCATE (XI, Yl) 

Yl-ABS(Yl) 

I F ( Yl.LT.3.0) Yl-0.0 
CALL NEWLIN2 

C ** GETS ROOT TRAILING POINT. 

CALL AOUTS T ( 20, A1 ) 

CALL AOUTST ( 19, TOP ) 

CALL RONEPNT ( X 2 , Y2 ) 

CALL LOC A T E ( X2, Y2 ) 

Y2«ABS(Y2 ) 

C0RD1-SQRT ( (Xl-X2)+*2+ ( Y1-Y2 )**?) 

C ♦* GETS VERTICAL LOCATION OF ROOT CHORD. 
CALL FRNTL0C(Z1,FTP,29 ) 
IF(YLOC.EQ.O.O)GO TO 12 
IF(IPAN3.LT.0)G0 TO 12 
Yl-WY (3 ) 

Z1 - WZ (3 ) 

12 CALL NEWLIN2 
IANSR2=I ANS ER 

C ** GETS LEADING POINT FROM PANEL TT P . 
CALL AOUTST (20, A2) 
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CALL AOUTS T ( 19» TOP ) 

CALL RONEPNT ( X3 * Y3 ) 

CALL LOCATE (X3> Y3) 

Y 3* ABS ( Y 3 ) 

IF(Y3.LT.3.0) Y3-0.0 
CALL NEWLIN2 

C ** GET TRAILING POINT FROM PANEL TIP. 

CALL AOUTST (19,A1B ) 

CALL AOUTST ( 19# TOP ) 

CALL RONEPNT ( X4> YA ) 

CALL LOCATE(XA, YA) 

YA* ABS ( Y A ) 

CORD2»SQRT( (X3-XA>**2+( Y3-YA>**2) 

C +* GETS VERTICAL LOCATION OF PANEL TIP. 

CALL FRNTLOC(Z2>FTP,29) 

CALL NEWLIN2 
CALL ANMODE 

C ** VARIOUS PARAMETERS ARE TO BE INPUT BY THE USER. 

PRINT "NUMBER OF SPANWISE VORTICES (NVOR)» 

1000 READ ( 1> * )N VQR 

IF (EOF (1 ) )1000> 1001 

1001 PRINT "NUMBER OF STREAMWISE VORTTCES (RNCV)" 

1002 READ ( 1# * ) RNC V 

I F ( EOF ( 1 ) ) 1002*987 

987 PRINT -^"LEADING EDGE SUCTION MULTIPLIER (SPC) 0 . 0<*S PC <« 1 . 0, 
PRINT *,"SPC<0.0 FOR POLHAMUS ANALOGY" 

1003 READ <1,* )SPC 

IF ( EOF < 1 ) > 1003> 100A 
100A IF (SPC .GT .1 .0 )G0 TO 987 
IF( IANSR2.EQ.87)GQ TO 6 

PRINT "ENTER THE ANGLE OF THE ROOT CHORO LINE (IN DEG.)" 
PRINT *>" A POSITIVE ANGLE CAUSES A DEFLECTION DOWN." 

1005 READ <1,*)AINC1 

IF ( EOF ( 1 ) )1005>1006 

1006 AINC1-TAN ( - A I NC 1 / 1 8 0 . * P I E ) 
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6 PRINT *> M ENTER ThF~ANGLE OF THF TIP CHORD LINE (IN DEG.)" 

1007 READ (1#*)AINC2 

IF ( EOF ( 1 ) ) 1007# 1008 

1008 AINC2-TAN ( - A I NC 2 /1 8 0 . * P I E ) 

I F < ITESTW.EQ.-l >G0 TO 78 

C ** WING CAMBER CAN BE ADDED. 

IF (NAP.EQ.O )G0 TO 7 
PRINT *#" " 

PRINT *#" " 

PRINT *#" " 

7 CALL WNGCAMB(NAP,ITESTW#CAMBl#CAMBO) 

Z1«Z1-C0RD1*CAMBI 

Z2«Z2-C0RD2*CAMB0 

78 IF(IPAN2.EQ.2)G0 TC 20 
I P AN2 *2 
WX(1)«X1 
WY(l)-ABSCYl) 

WZ(1)«Z1 

20 CONTINUE 
CALL NEWPAG 

C +* THE PANEL CAN BE ABORTED. 

PRINT * » "DO YOU WANT THE PREVIOUS PANEL" 

CALL Y ESNOF ( N02 ) 

IF(N02.LT.0)NAP«0 
IF (N02)87#21#21 

C ** VARIABLES TO SAVE DATA FOR ERROR RECOVERY • 

21 A X ■ X 1 
AY-Y1 
AZ-Z1 

ANG1-AINC2 

ACD*C0RD1 

AX2-SAFEX1-X3 

AY2-SAFEY1-Y3 

AZ2 *SAF EZ1 * Z2 

ACD2«SAFCRD1»C0RD2 
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ISAFEW-ITESTW 

SAFENAP*NAP 

CALL RESCALE (AXf AY >AZ>ACD> AX H»AY?»AZ?tACn?>0.0#l) 

C ** WRITES THE PANEL DATA UN A DISC FILE. 

WRITE (3> 100 ) AX# AY> AZ> ACD 
WRITE ( It 10G)AX2>AY2>AZ2#ACD2 
WRITE(3#101 )NVOR>RNCV>SPC>PDL 

WRITE(3>103)AINC1>AINC2>ITS>NAP>TQUANT*TSYNT>NPP 

NPANS-NPANS+1 
IPANEL-IPANEL+1 
I P ANA «-l 
IPAN3«-1 

IF (NAP.LE .2 )G0 TO 919 
WRITE(3#102)(XAF(J )>J*1>NAP) 

WRITE (3> 102 ) <ZC(1> J)> J«1*NAP> 
WRITE(3>102)(ZC(2>J)>J»1>NAP) 

DO 30 JJ“1»NAP 

JJ2-JJ+100 

XAF( JJ2) « X A F ( JJ ) 

30 ZC(2#JJ2)«ZC(2, JJ> 

919 CONTINUE 

100 F0RMAT(AF10.3»10X, "FLAT PANEL") 

101 F0RMAT( AF1C.C* 10X> »NVQR> RNC V> SPC > PDL " ) 

102 FORMAT ( 8F 10 • A ) 

103 F0RMAT(2F10.A>I2>8X>12>3(9X*I1)) 

GO TO 86 

87 IF < ISLT2 )86,85#85 
86 ISLT2 — 1 

C +* ALLOWS THE USER TO CHOOSE WHAT TYPE OF FLAT 
C *♦ PANEL IS DESIRED IF ANY. 

CALL NEWPAG 
CALL ANMODE 
NLINE-2 

IF ( .NOT • F IRST )G0 TO 2 

PRINT * t "ONE OF THREE POSSIBILITIES MUST BE CHOSEN NEXT*" 
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APPENDIX B. - Continued 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


**"( 1 ) 


**" 

*,"[2 ) 

*» " ( 3 ) 
* , n n 


A PANEL CAN BE ADDED 
PREVIOUS ONE (I.E. A 
AN ENTIRELY NEW FLAT PANEL 
THE NEXT TYPE OF PANEL CAN 


TO THE OUTER EDGE OF THE" 
MULTT-PANELFD WING)*" 

CAN BE ENTERED*" 

RE ENTERED." 


2 IF(FIRST)NLINE«8 
FIRST-. FALSE. 

PRINT **"FOR AN ENTENDING FLAT PANEL INPUT. ..E 

PRINT **"FQR A NEW FLAT PANEL INPUT N 

PRINT **"TG STOP FLAT PANELS INPUT ...0 

222 READ ( 1* 10A ) IANS 
IF (EOF (1 ) )2 22*223 

223 CALL L INF E E D ( NL INE ) 

IANSER-IANS(l) 

IF(IANSER.EQ.1HN)G0 
IF( IANSER.EQ.1HQJGC 
IF(IANSER.EQ.1HE)G0 


TO 
TO 
TO 

PRINT ** "RE-ENTER YOUR 


10 

81 

335 

ANSWER 


LEFT JUSTIFIED" 


(EXTENDING)." 
(NEW PANEL ) .» 
(QUIT) .» 


NLINE-2 
GO TO 222 

C ** HAKES IT POSSIBLE TG HAVE JOINED WING PANELS AND 
C ** RECOVER FROM ERRONEOUS INPUT. 

335 X1-SAFEX1 
Y1-SAFEY1 
Z1-SAFEZ1 
AINC1-ANG1 
C0RD1-SAFCRD1 
ITESTW-1SAFEW 
NAP-SAFENAP 
IF (NQ2.GT .0 )GQ TO 12 
DO 31 JJ-1*NAP 
JJ2-JJ+100 
XAF(JJ)-XAF(JJ2) 

31 ZC(2*JJ)«ZC(2»JJ2) 

GO TO 12 
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A 00 
A01 
A 02 
A 03 
AOA 
A05 
A06 
A07 
A 08 
A 09 
A10 
All 
A 12 
A 13 
A 1A 
A 15 
A 16 
A 17 
A 18 
A 19 
A 20 
A 21 
A 22 
A 2 3 
A2A 
A 25 
A26 
A 27 
A 28 
A29 
A 30 
A 31 
A 32 
A 33 
A 3A 
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81 

104 

C 

c ** 


c ** 


c ** 

40 


12 

C ** 

30 

C ** 


WX { 2 ) B S AF EX 1 

VORTAB 

435 

W Y ( 2 ) *A B S (SAFEY1) 

VORTAB 

436 

WZ (2 ) -SAFEZ1 

VORTAB 

437 

RETURN 

VORTAB 

438 

FORMAT ( 15 A1 ) 

VORTAB 

4 39 

END 

VORTAB 

440 


VORTAB 

441 

SUBROUTINE WNGC AMB ( NA P , I S L T> C AMB I , C A MBO ) 

VORTAB 

442 

ALLOWS WINGS TO HAVE CAMBER. 

VORTAB 

443 

COMMON / Z X P R / ZC(2>150)*XAF(150)>PHT(9),Rn(9) 

VORTAB 

444 

DIMENSION I ANS (15) 

VORTAB 

445 

IF (NAP.GT.C)GO TO 10 

VORTAB 

446 

CALL NEWPAG 

VORTAB 

447 

CALL ANMODE 

VORTAB 

448 

ASKS IF CAMBER IS DESIRED. 

VORTAB 

449 

PRINT *,"DC YOU WANT CAMBER IN THE PANELS” 

VORTAB 

450 

CALL YESNOF ( ISLT ) 

VORTAB 

451 

IF ( ISLT.GT.O)GO TO 40 

VORTAB 

452 

CAMBI-0.0 

VORTAB 

453 

CAMBO-O.O 

VORTAB 

454 

RETURN 

VORTAB 

455 

DETERMINES METHOD OF REPRESENTING CAMBER. 

VORTAB 

4 56 

PRINT *>”CAMBER IS ENTERED IN T-CHORD OR LINEAR MEASURE.” 

VORTAB 

457 

PRINT *,”D0 YOU WANT CAMBER IN T-CHORO" 

VORTAR 

458 

CALL YESNOF ( ICMB ) 

VORTAB 

459 

IF ( ICMB.GT.OJGO TO 12 

VORTAB 

4 60 

CALL WNGCMB2(NAP>CAMBI>CAMB0) 

VORTAB 

461 

RETURN 

VORTAB 

462 

PRINT " I NPUT CAMBER DATA BELOW.” 

VORTAB 

463 

ENTERS PANEL CAMBER IN 2-CHORD. 

VORTAB 

4 64 

PRINT ** "INPUT PERCENT CHORD < 0.0 TO STOP” 

VORTAB 

4 65 

NAP *NAP + 1 

VORTAB 

466 

ENTERS DATA FOR ROOT AND TIP OF PANEL. 

VORTAB 

467 


l f 


t * 
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PRINT *,"Z CHORD* ROUT CAMBER. E TTP CAMBFR" 

1009 READ ( 1* *') XAF (NAP ) >ZC ( 1* NAP ) * ZC ( 2*NA» ) 

IF (EOF (1 ) >1009*1010 

1010 IF (XAF (NAP ) .LT.O.OJGO TO 35 
GO TO 30 

35 NAP-NAP-1 

C ** CALCULATES PARAMETERS TO RELOCATF ROOT OF 
C ** PANEL IN VERTICAL DIRECTION. 

CAMBI-ZC (1*1) 

DO 20 J J ■ 1 » NAP 

20 IF ( ZC (1* J J ) .GT.CAMBI )CAMBI*ZC (1* J J ) 

CAMBI*CAMBI*.01 
GO TO 33 

C ** GETS CAMBER DATA IN LINEAR MEASURF. 

10 IF (ICMB.GT.O)GO TO 11 

CALL WNGCMB2(NAP*CAMBI*CAMBQ) 

RETURN 

C ** GETS DATA FOR PANEL TIP ONLY. 

11 PRINT *, "PERCENT CHORD AND ROOT CAMRER WILL BE DISPLAYED*" 
PRINT ** "FOR THE CURRENT PANEL." 

IANS ( 1 ) ■1000 
DO 150 J * 1 * NA P 
ZC(1»J)-ZC(2*J) 

WRITE 200*XAF (J >*ZC(1* J ) 

1011 READ ( 1* * ) ZC ( 2* J ) 

IF(E0F(1) >1011*150 
150 CONTINUE 

200 FORMAT ("INPUT TIP CAMBER ? Z CHORD *», F8 . 3 * 

1" ROOT CAMBER* "F7. 3*" Z") 

C ♦* RELOCATES PANEL TIP IN VERTICAL DIRECTION. 

C AMB1 *0 . 0 
33 CAMBO*ZC (2*1) 

DO 21 JJ«1*NAP 

21 IF(ZC(2*JJ>.GT.CAMBO)CAMBO«ZC(2»JJ> 

CAMB0«CAMBC*.01 
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A 88 
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A 89 
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VORTAB 

A 92 

VORTAB 
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A95 

VORTAB 

A96 

VORTAB 

A97 

VORTAB 
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A 99 
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502 
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tu TURN 
END 


VORTAB 

VORTAB 

VORTAB 


SUBROUTINE WNGCMB2 (NAP*CAMBI* CAMBO) 

C ** ALLOWS WINGS TO HAVE CAMBER* INPUT TN LINEAR DIMENSIONS. 
COMMON / Z X PR / ZC ( 2* 15C ) * X AF < 1 50 ) , PHT (9),P0(9) 

ITEST — 1 

IF ( NAP .EQ .0 ) GO TO 5 

ITEST-1 

GO TO 10 

C ** GETS PANEL ROOT DATA. 

5 PRINT "INPUT CAMBER DATA BELOW" 

PRINT * * " TO STOP INPUT DISTANCE FROM LEADING EDGE < 0.0" 
30 NA P * NA P + 1 

PRINT ** "ENTER DIST. FROM LEAD. EOGF » ANn ROOT CAMBER." 

1012 READ (1*+)XAF(NAP)*ZC(1*NAP) 

IF ( EOF ( 1 ) 11012*1013 

1013 IF (XAF (NAP ) .LT.0.0 )G0 TO 35 
GO TO 30 

35 NA P * NA P— 1 

C ♦* RELOCATES ROOT VERTICALLY. 

CAMBI«ZC( 1*1) 

DC 11 J J *2* NAP 

11 IF(ZC(1*JJ).GT.CAMBI)CAMBI»ZC(1*JJ) 

CAMBI*CAMBI/XAF (NAP) 

C ** PUTS ROOT DATA IN 7.-CH0RD. 

RC0RD1N*100./XAF(NAP) 

DO 12 JJ-1*NAP 
XAF(JJ)-XAF(JJ)*RC0RD1N 

12 ZC(1*JJ)-ZC(1*JJ)*RC0RD1N 
C ** GETS PANEL TIP DATA. 

10 DO 15 JJ-l.NAP 

PRINT ** "STATION NUMBER "»JJ 
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503 
5 OA 
505 


506 

507 

508 
5 09 
5 10 

511 

512 

513 
5 1 A 
5 15 
5 16 

517 

518 

519 
5 20 
521 
5 22 
523 
52A 
5 25 

526 

527 

528 
5 29 
5 30 
531 
5 32 
533 
53A 
5 35 
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PRINT *# "ENTER DIST. FROM LEAD. PDGF fi TIP CAMBER." 
101 A READ (1#*)XAF(JJ)#ZC(2#JJ) 

I F ( EOF ( 1 ) ) 101A# 15 

15 CONTINUE 

IF{ITEST.GT.O)CAMBI*O.C 
C *'* RELOCATES PANEL TIP VERTICALLY. 

33 CAMBO-ZC (2#1 ) 

C ♦♦ PUTS DATA IN t-CHORD. 

DO 16 JK«2#NAP 

16 IF ( ZC (2# JK ) .GT.CAMBO)CAMBO*ZC <?» JK ) 

CAMBO-CAMBO/XAF (NAP ) 

RC0RD2N«100./XAF(NAP ) 

DO 17 JK * 1 # NA P 

17 ZC (2# JK) *ZC (2, JK)*RC0RD2N 
RETURN 

END 


SUBROUTINE YESNOF ( WHE R E ) 

C ** ALLOWS A QUESTION TG BE ANSWERED WITH YES OR NO# 

C ** OR# TRUE OR FALSE# CR ANYTHING L p SS THFN 16 CHARACTERS 
C ** THAT STARTS WITH Y»T»N» OR F. 

INTEGER WHERE 
DIMENSION ANS ER ( 15 ) 

PRINT *# "ENTER Y OR T FOR YES# OR N OR F FOR NO" 

A READ <1#1) ANSER 
I F ( EOF ( 1 ) 

6 CALL LINEF 
GO TO A 

1 FOR MAT ( 1 5 A 1 ) 

10 ANS *ANSER ( 1 ) 

IF (ANS.EQ.IHT)GO TO 2 
IF(ANS.EQ.1HY)G0 TO 2 
IF (ANS.EQ.IHN)GO TO 3 
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IF (ANS.EQ.1HF )60 TO 3 

PRINT *# "RE-ENTER YOUR ANSWER' LEFT J IJSTTFT FD • " 
GO TO 4 
C ** YES 

2 WHERE«1 
RETURN 

C ** NO 

3 WHERE *— 1 
RETURN 
ENO 

C 


SUBROUTINE VTAIL 

C ** GENERATES VERTICAL FLAT PANELS. 

COMMON /FUSS/XO# YO# FSCORD# F S L OP F , BFUS#NPANS#FSCALE 
DIMENSION S IDEV (2 ) 

DIMENSION AO(3)#A2(2)#TOP<5),Aim#AlB(?) 

REAL NVOR 

DATA SPC# A I NCI# AINC2# ITS/ NAP, IQ1JANT# I S YN T , N PP / 3*0 . 0 . 5 *07 
DATA PDL/O.O/ 

DATA TOP / 10HA POINT F0#10HR SPANWISFrlOH LOCATION 
1 > 10H ( FROM THE #10HT0P VIEW) / 

DATA SIDEV/1CH(FR0M THE #10HSIOF VTFW)/ 

DATA A0/10HB0TT0M LEA# 10HDING PniMT,lH / 

DATA A1/10HBQTT0M TRA# 10HIL ING POIN,?MT / 

DATA A1B / 1CHT0P TRAIL I#9HNG POINT / 

DATA A2/10HT0P L EADIN# 10HG POINT / 

C ** ASKS IF VERTICAL PANELS ARE DESTPEO. 

CALL NEWPAG 
CALL ANMODE 

PRINT *#"D0 YOU WANT VERTICLE PANE! S'* 

CALL YESNOF(ISLT) 

IF( ISLT.LT.OJRETURN 
110 CALL LINFEtD(4) 
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16 


C ♦♦ 


c ** 


c +* 
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c ** 


APPENDIX B. - Continued 


GETS BOTTOM LEADING POINT. 

CALL AOUTS T (2 1, AO ) 

CALL AOUTST (20,SIDEV) 

CALL RONEPNT(Xl, Z1 ) 

CALL LOC A TE S ( XI, Z1 ) 

CALL NEWLIN2 

GETS THE BOTTOM TRAILING POINT. 

CALL AOU TS T ( 22, A1 ) 

CALL AOUTST (20, SIDEV) 

CALL RONE PNT ( X2, Z2 ) 

CALL LOC A TES ( X2, Z2 ) 

CORD1 *SQR T ( (X1-X2)**2+(Z1-Z2>**?> 

CALL NEWLIN2 
GETS TOP LEADING POINT. 

CALL A0UTST(18,A2) 

CALL AOUTST (20, SIDEV) 

CALL PON E PN T ( X3 , Z3 ) 

CALL LOC A TE S ( X3» Z3 ) 

CALL NEWL1N2 

GETS TOP TRAILING PC1NT. 

CALL A0UTST(19,A1B) 

CALL AOUTST (20, SIDEV) 

CALL RONE PN T ( X4 , Z4 ) 

CALL LOC AT ES ( X4 , Z4 ) 

CORD2 *SQRT ( ( X 3-X4 ) **2+ ( Z3-Z4 ) **2 ) 

CALL NEWLIN2 

LOCATES PANEL IN SPANWISE DIRECTION • 

CALL AOUTST (50, TOP ) 

CALL RQNEPNT(AX,AY) 

CALL LOCATE (AX, AY) 

Yl-AY 

IF(ABS(Y1).LT.4.0)Y1*0.0 

GETS THE NUMBER OF SPANWISE AND CHOP DMT S F VORTICES. 
CALL NEWLIN2 
CALL ANMODE 
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PRINT ** "NUMBER OF SPANWISE VOOTTCCS (NVOR)" 

VORTAB 

637 

1016 

READ ( 1 » ♦ ) N VQR 

VORTAB 

6 38 


IF(EOKl) ) 1016* 1017 

VORTAB 

639 

1017 

PRINT ** "NUMBER OF STREAMW1S E VORTICES (BNCV)" 

VORTAB 

6 AO 

1018 

READ ( 1 * * ) R NC V 

VORTAB 

641 


IF ( EOF ( 1 ) ) 1018* 1019 

VORTAB 

642 

1019 

CALL LINFEED(6) 

VORTAB 

643 

c ♦♦ 

ALLOWS THE PANEL TO BE ABORTED. 

VORTAB 

644 


PRINT * * " D 0 YOU WANT THE PREVIOUS PANEL" 

VORTAB 

6 45 


CALL YfcSNOF (NO ) 

VORTAB 

6 46 


I F ( NO ) 11# 2 1 * 2 1 
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647 

21 

AX-X1 
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AX2«X3 
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AY-AY2-Y1 
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AZ-Z1 
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AZ2 «Z3 
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ACD-C0RD1 
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ACD2-C0PD2 
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654 


CALL RESCAL E ( AX* A Y* AZ * AC D> AX2 »AY? , A 7 2 > AC 07 * 0 .0 * 1 ) 

VORTAB 

655 

C *♦ 

WRITES THE PANEL DATA ON DISC. 

VORTAB 

656 


WRITE(3*100)AX*AY* AZ* ACD 
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WRITE(3*100)AX2#AY2#AZ2»ACD2 
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6 58 


WRITE (3*101 )NVOR*RNC V*SPC> POL 

VORTAB 

6 59 


WRITE (3* 102 IAINC1* AINC2* I TS * NAP , I QUANT * I S YNT » NPP 

VORTAB 

660 


NPANS “NPANS +1 

VORTAB 

661 

100 

FORMAT(4F10.3*10X* "VERTICAL PANEL") 

VORTAB 

662 

101 

FORMAT (4F10.Q* 10X, "NVOR*RNCV*SPC»PDL") 

VORTAB 

663 

102 

FORMAT (2F 10.2* 12*8 X* 12* 3 (9X, ID ) 
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11 

CALL NEWLIN2 

VORTAB 

665 

C ** 

ALLOWS THE USER TO ADD ANOTHER PANFL I P DESIRED. 

VORTAB 

666 


PRINT *>"D0 YOU WANT ANOTHER VERTICAL PANEL" 

VORTAB 

667 


CALL YESNOF(ISLT) 

VORTAB 

668 


IF(ISLT.LT.O)RETURN 

VORTAB 

669 


CALL SYSTEM(0*" ") 

VORTAB 

6 70 


CALL NEWPAG 

VORTAB 
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GO TO 16 


— 

VORTAB 

672 


END 
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c 
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SUBROUTINE NACELLE 

VORTAB 

675 

C 

t* 

CALCULATES DATA FOR NACELLE PANELS. 

VORTAB 

6 76 



COMMON /FUSS/XO, YO*FSCORD»FSLCPF, BFUS,NP ANS>FSCALE 

VORTAB 

677 



COMMON / Z X P R / ZC(2,15C),XAF<150>»PHT(9),P0(9) 
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REAL NVOR 

VORTAB 

6 79 



DIMENSION FTP(2),PHI2<9) 

VORTAB 

6 80 



DATA SPC,AINC1> A INC 2# I TS > I QU ANT ♦ I S YNT , NP P / 3*0 . 0, 3*0, 1 / 

VORTAB 

681 



DATA PHI 2 /-90.0, -4 6. 0,0.0 , 45. 0,90.0,135 .0, 180.0,225.0,270.0/ 

VORTAB 

682 



DATA FTP/1 CHTOP NACELL, 8HE POINT / 

VORTAB 

683 

C 

** 

REINITIALIZES PHI ARAY. 

VORTAB 

684 



DO 1 J*l,9 

VORTAB 

685 


1 

PHI ( J)*PHI2 < J ) 

VORTAB 

686 

C 


ASKS IF NACELLES ARE DESIRED. 

VORTAB 

687 



NLINZ-5 

VORTAB 

6 88 



CALL NEWPAG 

VORTAB 

689 



CALL ANMODE 

VORTAB 

690 



PRINT *,»'DG YOU WANT NACELLES" 

VORTAB 

691 



CALL YESNDF ( ISL T ) 

VORTAB 

692 



IF(ISLT.LT.O)RETURN 

VORTAB 

693 

C 

** 

FRONT TWO NACELLE POINTS. 

VORTAB 

6 94 


5C 

1 PRINT 104 

VORTAB 

695 



CALL SYSTEM (0 * " ") 
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696 



CALL LINFEED(NLINZ ) 
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697 



CALL RONE PNT ( X 1# Y1 ) 
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698 



CALL LOCATE (XI, Yl) 
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6 99 



PRINT 105 

VORTAB 

700 



CALL SYST EM (0, " ") 
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701 



CALL NEWLIN2 

VORTAB 

702 



CALL RUNEPNT(X,Y) 

VORTAR 

703 



CALL LOCATE (X,Y) 

VORTAB 

704 
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R«ABS(Y1-Y)*.5 

YCENTER«(Y+Y1)*.5 

IF(ABS(YCENTER).LT.5.)YCENTER-0.0 

NAP-1 

ZC ( 2 # 1 ) -R *R 
XAF(1)=ZC <l,l>-0.0 
CALL NEWLIN2 

C ** INSTRUCTIONS ARE DISPLAYED. 

CALL ANMODE 

PRINT "ENTER BODY POINTS ALONG ONE SIDF ( E XC LUD I NG FRONT " 
PRINT "POINTS). PROGRESSING FROM FRONT TO REAR AFTER THE" 
PRINT *>"LAST POINT HAS BEEN PROCESSED FNTER Q FOLLOWED BY" 
PRINT +> "A CARRIAGE RETURN." . 

CALL S YST EM (0> " ") 

CALL LINFEED(3) 

C *♦ BODY POINTS READ FROM TABLET. 

9 CALL TAB ARM ( 0 ) 

CALL GE TP NT ( I » I X > I Y ) 

X-FLOAT ( I X ) 

Y- FLOAT ( I Y) 

I F ( I • EQ • 8 1 ) GO TO 11 

IF ( I. GT. 31) GO TO 9 

IF(X.EQ.O.O.A.Y.EQ.O.O)GO TO 9 

CALL LOCATE C X > Y 

NAP-NAP +1 

R2-Y-YCENTER 

ZC(2>NAP)-(R2*R2) 

XAF(NAP) -X-X1 
GO TO 9 

11 CALL TABOFF 
C0RD1-XAF (NAP) 

RCEP1-100.C/C0RD1 
DO 14 J »2>NAP 
ZC(1> J)-0.0 

1 A XAF(J )«XAF( J)*RCEP1 


X 4 
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705 

706 

707 

708 
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7 1 A 
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716 
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PDL-999. 

c ** TOP OF NACELLE IS LCCATED VERTICALLY. 

CALL FRNTLCC(Z»FTP,18> 

CALL NEWLIN 
CALL ANMODE 
CALL SHI FT ( NAP ) 

Z2«Z-2.0*R0( 1 ) 

NV0R-6.0 

IF(YCENTER)5,4,5 
A NV0R«4 .0 
& NV0R1-INT(NV0R)+1 
IF(NV0R.EQ.8)Z«Z2 

C ** NUMBER OF STREAMWISE VORTICES IS DETERMINED. 

PRINT 106 

30 READ( 1,+ IRNCV 
IF ( EOF ( 1 ) >30,31 

C ** ALLOWS PANELS TO BE ABORTED. 

31 PRINT *,"DG YOU WANT THE PREVIOUS PANEL" 

CALL YESNOF(NO) 

IF( NO >99,20,20 
20 NPANS*NPANS+1 
AX«AX2«X1 
AY»AY2“YC ENTER 
AZ»Z2 
AZ2«Z 

ACO-ACD2-CORD1 

CALL RESCALE(AX,AY,AZ,ACD,AX2,AY2,A72, AC02,999.0,NV0R1) 
C ** WRITES PANEL DATA ON DISC. 

WRITE (3,100>AX,AY, AZ, ACD 
WRITE (3, 100>AX2, AY2,AZ2,ACD2 
WRITE (3,101 )NVUR,RNCV,SPC,POL 
WRITE (3, 102 > (PHI ( J > ,R 0 ( J > , J* 1 >NV0P1 ) 

WRITE (3, 103) A INC 1, A INC 2, ITS, NAP, TOIJANT » T S YNT » NPP 

IF(NAP.LE .2 >00 TO 99 

WRITE (3,102 > (XAF(J ),J-1, NAP ) 
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Ol 

ro 

WRITE (3, 102)<ZC(I, J),J-1,NAP) 

WRITE (3, 102 ) ( ZC (2, J ), J-1,NAP ) 

C ** ALLOWS ADDITIONAL NACELLE PANELS TO 8F FNTFRFD. 

99 PRINT *,»IS ANOTHER NACELLE DESTRFD" 

CALL YESNOF ( ISL T ) 

IFtlSLT .LE.O)RETURN 
CALL NEWPAG ' 

NL INZ-2 
CALL ANMODE 
GO TO 50 

100 FORMAT (4F10.3,10X, "NACELLE PANEL") 

101 F0RMAT(4F10.0,10X, "NVOR,RNCV,SPC,pn|.") 

102 F0RMAT(8F1C.4) 

103 F0RMAT(2F10.2,12,8X,I2,3(9X, ID ) 

104 FORMAT ( " FROM THE TCP VIEW'S/A 

1" INNER OR OUTER FRONT NACELLE POINT") 

105 FORMAT!//," REMAINING FRONT NACFLLF POINT") 

106 FORMAT (/, "NUMBER OF CHORDWlSE VORTTCES (RNCV)») 
END 
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VORTAB 781 
VORTAB 782 
VORTAB 783 
VORTAB 784 
VORTAB 785 
VORTAB 786 
VORTAB 787 
VORTAB 788 
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SUBROUTINE SHIFT(NAP) 

C ** LOCATES NACELLES IN THE VERTICAL DIRECTION. 

COMMON / Z X PR / ZC (2, 150) ,XAF ( 160) , PHT (9) ,R0 (9) 
TE ST-ZC ( 2 , 1 ) 

DO 10 1-2, NAP 

10 IF(ZC(2,1).GT.TEST)TEST=ZC(2,I) 

R-SQRT (TEST ) 

DO 11 J * 1 > 9 

11 R0( J ) *R 

- REFARE A- TEST*. 01 
DO 12 J ■ 1 , N AP 

12 ZC(2, J)-ZC(2, J)/REFAREA 
RETURN 
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808 
8 09 


C 
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SUBROUTINE FU S L AGE ( YL OC ) 

C +* THIS GENERATES THE FUSELAGE PANFLS. 

COMMON / FUS S / XO, YO,FS CORD, FS LOPE » «FIIS,N PANS, FSC ALE 
COMMON /FUSS2/RN0T,NAP 

COMMON /ZXPR/ ZC(2,150),XAF(150)»PHT(9)»R0(9) 

COMMON / Z S E A / ZC2 ( 2 > 5 0 ) , X A F 2 ( SO ) 

REAL NVOR , KORD 
DIMENSION PHI 2 ( 9 ) 

DATA SPC,AINC1,AINC2, ITS# I QU ANT » T S YNT ♦ NP P / 3 *0 . 0, -1 » 0, 0, 1 / 

DATA PHI2/-90. 0,-4 5. 0,0. 0*45. 0,90. 0*135. 0,1 80. 0,225. 0,270.0/ 

C +* GETS NOSE POINT FROM TOP VIEW. 

50 PRINT *, "ENTER NOSE POINT FROM TOP VIFW" 

CALL SYSTEM(0»" ") 

CALL LINFEED(4) 

CALL RGNEPNT! ANYX, YLOC > 

CALL LOCATE (ANYX, YLOC ) 

IF(YLC)C.LT.3.0)YL0C»0.0 
YLOC* ABS ( YLOC ) 

CALL NEWPAG 
NO * 1 

CALL L1NFEEDI2) 

PDL2 *999 . 

C ** INITIALIZES PHI ARRAY. 

DO 1 J • 1 , 9 
1 PHI ( J)*PHI2 ( J) 

C +* INSTRUCTIONS ARE TO BE PRINTED. 

NLINES-6 
CALL ANMODE 
' PRINT 321 

321 FORMAT!" TO MAKE A FUSELAGE BODY ENTER POINTS FROM THE SIDE"* 

1" VIEW. THE NOSE POINT IS FIRST. THE NEXT POINTS APE " 
2"ENTERED IN",/," PAIRS, THE FIRST FROM THE TOP OF THE FUSELAGE" 
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8 41 
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3"# THE SECOND FROM THE BOTTOM DIRECTLY BHOW THE FIRST ."# /# 
4" THE LAST TWG COME FROM THE END OF THE FUSELAGE# AFTER'* 

5" THESE TWO ENTER 0 TO END"#//#" NOSE POTNT") 

C ** GETS NOSE POINT. 

41 CALL SYSTEM(0»" "> 

CALL LINFEED(NLINES ) 

CALL RONE PNT ( XI# Z1 ) 

CALL LOCATES ( XI# Z1 ) 

IF (ABS(Z1 ) .LT.3.0)Z1«C.0 
ZLOC-Z1 

IF ( ABS (X1).LT.3.C)X1»C.G 

XAF1-X1 

XAF(l)-0.0 

ZCd#l)«0.0 

ZC ( 2# 1) *0 ♦ 0 

NA P • 1 

C *+ GETS BODY POINTS. 

90 PRINT *#" TOP POINT" 

CALL SYSTEM(0#" ") 

CALL NEWLIN 
9 CALL TABARM (0 ) 

CALL GETPNT ( I# IX# I Z ) 

X»FLOAT( IX ) 

Z-FLOAT(IZ) 

IF(I.E0.81)G0 TO 11 
IFd.GT.3DG0 TO 9 
CALL LOCATES(X#Z) 

PRINT *#" BOTTOM POINT" 

CALL SYSTEMtO#" "> 

CALL NEWLIN 
209 CALL TABARM(O) 

CALL GETPNT ( 1 # I X 2# IZ?) 

IF(I.EQ.81)GU TU 11 
IF(I.GT.31)GG TO 2C9 
X2«FL0AT( 1X2) 

Z2-FL0AT(IZ2) 
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VOR TAB 

847 
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VORTAB 
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R 5 3 
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854 
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CALL LOCATES (X2#Z2 ) 

NAP-NAP+ 1 

XAF (NAP) -X-XAF1 

ZC (2>NAP ) « ( (Z -Z2 ) + .5 ) *+2 

ZC(1»NAP)*(Z + Z2 )* • 5— ZLQC 

GO TO 90 

11 CALL TABOFF 
CALL NEWLIN2 
CALL ANMDDE 

C ** ALLOWS DATA TO BE RE-ENTERED IF AN ERROR HAS BEEN MADE. 
PRINT *>"D(J YOU WANT TO RE-ENTER THIS DATA” 

CALL YESNOF (LBODY) 

IF ( LBODY .LT .0 )G0 TO AO 

PRINT "RE-ENTER FUSELAGE DATA." 

CALL NEWLIN 
GO TO bO 

AO I2-INT ( .5+NAP ) + l 

C ** DATA IS PUT IN TERMS CF %-CHURD. 

C0RD1-XAF (NAP ) 

CRDNRM-100.0/CURD1 
RNOT*SQRT (ZC(2,I2) ) 

RN0RM«100.C/ZC(2*I2) 

DO 12 J«1>NAP 
ZC(1»J)*ZC(1>J)*CRDNRM 
XAF(J)«XAF( J)*CRDNPM 

12 ZC(2/J)*ZC(2»J)*RNCRM 

C ** RADII OF FUSELAGE ARE DETERMINED. 

DO 13 J B 1 > 9 

13 R 0 ( J ) *RNOT 
PDL-999. 

Z2 *Z = -RNOT 
NVOR-8. 0 
IF(YL0C)5»A,‘j 

A NVOR-A 

Z«Z2+2.0*RN0T 
5 NV0R1-1+1NT (NVOR ) 
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8 78 

879 

880 
8 81 
882 
883 
8 8A 

885 

886 

887 

888 
8 89 

8 90 

891 

892 

893 
89A 

895 

896 

897 

898 

899 

900 

901 
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903 
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906 

907 

908 

909 

910 

911 

912 

913 
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IF(NVOR.EQ.8)Z-Z2 

** FUSELAGE CAMBER CAN BE ADDED HEPF. 

CALL RELOCATt YLOC» ZLOC» CORD1 , XI ) 

Z2*Z2+ZLQC 

Z-Z+ZLOC 

** ASKS FOR NUMBER OF FUSELAGE PANFLS. 

55 PRINT "NUMBER OF FUSELAGE PANELS" 

56 RE AD ( 1 > * ) NUM 

IF (EOF ( 1 ) )58*57 
58 CALL LINEF 
GO TO 56 

57 IF (NUM.GT.O.A.NUM.LT.5)GO TO 10 

PRINT "RE-ENTER YOUR NUMBER RFTW CC N 0 AND 5" 
CALL LINE EED ( 3 ) 

GO TO 55 

** ALLOWS MULTIPLE FUSELAGE PANELS. 

10 IF (N0.EQ.1)NPANS*NPANS+NUM 
XAF2 ( 1 ) *X AF ( 1 ) 

ZC2(1*1)>ZC (1*1) 

ZC2(2*1)«ZC<2>1) 

ANOT-O.O 

N0*1 

JTWO-1 

NLINZZ-7 

** LOOP TO MAKE MULTIPLE FUSELAGE PANFLS. 

DO 30 I MP ■ 1 #NUM 
54 JNOT *JTWO 

IF( IMP.EO.NUM)GO TO 22 
PRINT 10 A * I MP 

104 FORMATl" ENTER BREAK POINT FOR p ANFL # ">T1* 

1" » (FROM THE TCP VIEW)") 

CALL SYSTEMtO," ") 

CALL LINFEEDtNLINZZ) 

NLINZZ«8 

CALL RONEPNTl AX* AZ ) 

CALL LOCATESt AX,AZ) 
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914 

915 

916 

917 
9 18 

919 

920 

921 

922 

923 

924 

925 

926 

927 

928 
9 29 

930 

931 

932 

933 

934 

935 

936 

937 

938 
9 39 

940 

941 

942 

943 

944 

945 

946 

947 

948 

949 
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AX-AX-XAFl 
GC1 TO 23 

22 AX*XAF (NAP )*C0RD1*.01 
J T WO*NAP 

23 KORD*AX-ANOT 

IF( IMP.EQ.l )GO TO 36 

C ** FIRST POINT OF PANEL IS DETERMINED. 

ZC2( lfl) *ZC2( 1>NAP1 )*ULDCORD/KOPn 
Z C2(2,1)-ZC2(2,NAP1) 

X AF 2 ( 1 ) * 0 • 0 
36 Xll-ANOT 

OLDCOPD-KORD 
SAF ANOT 3 ANOT 
S AF J T WO* J T W 0 
ANOT-AX 

C ** NUMBER OF CHORDWISE VORT1CFS I S D FT FP MI NFD . 

IF ( IMP ,NE .NUM ) PR INT 106 
106 FORMAT (/ ) 

PRINT *, ‘'NUMBER UF CHORDWISE VHPTICTES (RNCV)" 

PRINT * > " F 0 P PANEL NUMBER ">LMP 

28 READ ( 1> * ) PNCV 
IF ( EOF ( 1 ) )60» 29 

60 CALL LINEF 
GO TO 28 

29 PCORD-lOO .+AX/C0RD1 

IF (IMP.EQ.NUM)GO TC 26 
00 31 J 3 J NO T * NAP 

31 IF(XAF( J) .GT.PCORDJGO TO 32 

32 JTWO* J 

C ♦* DATA AT PANEL JOINING POINTS ARF CALCULATED • 

NAPP-JTWO-JNOT+2 

IF(IMP.EQ.1)NAPP 3 JTW0 

RATIO*(PCGPO - XAF(J-1) )/ ( XAF ( J)-XAF(J-l) ) 

ZC2 (1*NAPP) -RATION (ZC ( 1» J )-ZC (1 »J-1 ) ) + ZC ) 

ZC2(2»NAPP)“RATI0*(S0RT (ZC(2»J))-SORT(7C(?»J-l )) ) +SORT ( ZC C 2 1 J.-l ) ) 
ZC2(2»NAPP)*ZC2(2»NAPP)+*2 
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XAF2 (NAPP ) >PCOKD 
GO TC 27 

++ CALCULATES THE PARAMETERS TO MUDTFY DATA FDR THF 
** APPROPRIATE PANELS. 

26 NAPP*JTWO-JNOT+2 

IF (NUM.EQ .1 )NAPP*NAP 

27 ZC 1FUG*C0RD1/K0RD 
XAFFUG1*ZC1FUG 
XAFFUG2--100.+X11/K0RD 
IF ( IMP . EO . 1 )X AFFUG2-0 .0 
IF ( IMP.t 0 . NUM ) GO TO 37 
ZC2(1>NAPP)»ZC2(1>NAPP)*ZC1FUG 

XAF2 (NAPP J-XAF2 (NAPP) ♦ X AF FUG 1 +X AF FUG? 

37 JNO* JNOT-1 
J 1 W * J T WO- 1 
NAP1-NAPP 

IF ( IMP. EQ.NUM ) JTW= JTWCi 
IF( IMP.GT.l ) JNOJNC-1 

IF(NAPP.EQ.2.A.IMP.LT.NUM)G0 TO 333 
*+ MODIFIES DATA FOR THE A P PROP R I ATF p A M F L AND WRITES IT 

** ON DISC. 

DO 33 JINC*JNOT» JTW 
J FUG ■ J INC - J NO 

ZC2(I>JFUG)°ZC(1>J1NC)*ZC1FUG 
ZC2 12 * JFUG ) «ZC ( 2 * J INC ) 

33 XAF2 (JFUG) «XAF( J INC )*XAFFUG1 +XAFFUG? 

333 CONTINUE 

** ALLOWS A PANEL TO Bt ABORTED. 

PRINT *>"D0 YOU WANT THE PREVIOUS PANF|.» 

CALL YESNOF (NO) 


IF 

{ NO ) 55 2 > 2 0> 

20 


IF 

(NAPI.F0.2) 

CALL NAP 

E 2 (NAP1 ) 

VAR 

IABLES APE 

PREPARED 

FOR RE-SCALING. 

AX 

•AX2-X1 1+XA 

FI 


AY 

-AY2-YL0C 



AZ 

= Z2 
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A/ 2 ■ Z 

VDRTAB 

1022 


AGO" ACD2 "KURD 

VDRTAR 

1023 


CALL RESCALE (AX,AY,AZ>ACD, AX2,AY7, A7’, ACD2 » PDL2,NV0R1 ) 

VORTAB 

102A 


PDL2*O.0 

VOP TAB 

1025 

** 

WRITES PANfcL DATA ON DISC 

VORTAB 

1026 


WRITE <3> 100)AX>AY, AZ>ACD 

VORTAB 

1027 


WRITE(3> 100)AX2*AY2»AZ2>ACD2 

VORTAB 

102B 


WRITEI3, 101 )NVOR,RNCV>SPC»PUL 

VORTAB 

1029 


WR I TE ( 3» 102 ) (PHI (J ),RG( J ) » J»1,NVCIP1 ) 

VORTAB 

1030 


WRITE (3# 103)AINC1» AINC2# ITS, NAP] » TOMANT» IS YNT,NPP 

VORTAB 

1031 


WRI TE (3> 102) <XAF2( J )* J-1,NAP1 ) 

VORTAB 

1032 


WRITE (3, 102) (ZC2<1, J)» J-1/NAP1) 

VORTAB 

1033 


WRITE ( 3> 102 ) ( Z C 2 ( 2 > J ) > J * 1> NAP 1 ) 

VORTAB 

103A 


GO TO 30 

VORTAB 

1035 

552 

JTWO-SAFJTWO 

VORTAB 

1036 


ZC 2 ( 1> NA P 1 )=ZC2 ( 1» 1 ) 

VORTAB 

1037 


ZC2(2,NAP1 )«ZC2(2,1) 

VORTAB 

1038 


ANOT-SAFANOT 

VORTAB 

1039 


GO TO 5A 

VORTAB 

10 AO 

30 

CONTINUE 

VORTAB 

10A1 

100 

F0RMAT(AF10. 3» 10X, "FUSELAGE PANEL") 

VORTAB 

10A2 

101 

FORMAT (<» F10. 3 ♦ 10X, "NVOR, R NC V, SPC » P DL " ) 

VORTAB 

10 A3 

102 

FORMAT (8 F 10 . A ) 

VORTAB 

10 AA 

103 

FORMAT (?F 10 • 2 » I 2 *8 X» 12* 3 ( 9X * ID ) 

VORTAB 

10A5 


IF(YLOC.NE.O)RETURN 

VORTAB 

10A6 


PRINT *,"IS ANUTHER FUSELAGE RFOIJIPPD” 

VORTAR 

10A7 


CALL YESNOF(INOTHER) 

VORTAB 

1 0 AH 


CALL S YS T EM ( 0» " ") 

VORTAB 

10A9 


CALL NEWPAG 

VORTAB 

1050 


CALL ANMODE 

VORTAB 

1051 


IF ( INUTHER .LT.O JRETURN 

VORTAB 

1052 


NO * 1 

VORTAB 

1053 


GO TO 50 

VORTAB 

10 5 A 


END 

VORTAB 
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1056 
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SUBROUTINE NAPE2(NAP1) 

C ++ CHANGES UNt SEGMENT PANELS INTO TWO SFGMFNT PANFLS . 
COMMON / Z S E A / ZC2 < 2>50 ) , XAF2 ( 50 ) 

NAP 1 ■ 3 

XAF2 (2)-50.0 
XAF2 (3) *1CC.C 
ZC2(1>3)»ZC2(1>2) 

ZC 2 ( 2# 3 ) 1 ZC 2 ( 2> 2 ) 

ZC2{1>2)«(ZC2{1#3)+ZC2(1>1)M‘.‘> 

S1*SQRT(ZC2 (2» 3 ) ) 

S2»SQRT<ZC2 (2,1 ) ) 

ZC2<2#2)*{ <Sl+S2)*«5)*+2 
RETURN 
END 
C 


SUBROUTINE R ELQC AT ( YL CC » ZLOC » FS C ORO » XS T A RT ) 

C ** ALLOWS WINGS TO BE ATTACHED TO THF FUSELAGE. 

COMMUN / Z XPR / ZC(2>150)>XAF(150)»PHT (9)tR0(9) 
COMMON /FUSS2/RNQT > N AP 

COMMON /WINS/ AWX(3),AWY(3),AWM3I.MT) = 
JIMENSIQN TES T ( 3 ) 

DATA FIE/3.13159/ 

IF(M0DE.EQ.-1)RETURN 

C ** LOCATES THE STARTING POINT OF THF WING WITH 
C ** RESPECT TO THE FUSELAGE. 

IWWW«1 

I F ( YLOC . E 0 . 0 ) GO TO 1 
IWWW*2 

1 WX«AWX(IWWW)-X START 
WY-AWY(IWWW) 

WZ * A WZ ( 1WWW ) 

AWZ ( 3 ) *W Z 

PC0RD«(100.O*WX)/FSCCRD 
WY-ABS (WY-YLOC) 
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IF ( YLOC .GT .0) AWY (3 ) -YLUC+WY 

C ** DETERMINES WHICH X-CHORD STATIONS STARTING POINT OF 
C ** WING IS BETWEEN. 

DO 2 J *2 > NA P 

2 IF (XAF ( J ) .GT.PCQRO )G0 TO 3 

3 RATIO(PCORD-XAFU-l) )/(XAF( J)-XAF( J-l) ) 

2SHIFT *ZC ( 1 > J-l ) + ( ZC( 1 * J ) -ZC ( 1> J“1 ) ) *R A T 1 0 
DO 10 IJK»1,NAP 

10 ZC(1,1JK)-ZC(1,IJK)-ZSHIFT 
ZL0C«ZLGC+ZSHIFT*.01*FSC0RD 
WZ-WZ-ZLOC 

RAD1-SQRT (ZC(2> J-l) )+ ( SORT < ZC (? * J) )-S0RT(ZC<2> J-l) ) )*RATIO 

RAD1 X RAD1*.1 

RADIUS*RAD1*RN0T 

C *+ DISTANCES BETWEEN WING 0 FUSELAGE SUBPANELS ARE DETERMINED. 

11 DO 30 I “ 1 > 3 
FI*Pit/180.0*PHl (1+1) 

ZTEM'SIN < FI )*RADIUS 
YTEM-COS ( FI )*RADIUS 

30 TESTt I ) «= S OR T ( < W Z-ZT EM ) **2 + ( W Y-YTEM ) **2 ) 

TESTI-TEST (1) 

TEST2-TEST (2 > 

TEST3*TEST(3) 

C *♦ SUBPANEL 60UN0RY CLOSEST TO THE WTNG IS SELECTED. 
IF(TESTl.LT.TEST2.A.TESTl.LT.TFST3)Gn TO AO 
IF(TEST2-TEST3)38>3B,39 

38 ITEST-3 
GO TO A 1 

39 ITEST-4 
GO TO A 1 

AO ITEST-2 

C ♦ * SHIFTS A SUBPANEL HOUNDKY TO CORRESPOND TO WING ROOT. 

A 1 R0TES1»SURTCWZ*WZ+WY*WY) 

R0( I TES T ) -ROTES T2»kOTES T/RAD1 

PHI(ITEST)-PHIT»ASIN(WZ/RL)TEST)*1B0.0/PIF 

DELTARO- (R0TEST2-RN0T ) 
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c 

** MAKES SIZE Oh SUBPANELS ARROUND FUSELAGE 

MORE EQUAL. 

VORTAR 

1127 


IF(ITtST-3>50, 51*52 


VORTAB 

1128 


50 PHIlNC-(90.0-PHIT) /3.0 


VORTAB 

1129 


DELT »DELT ARO/3 . 0 


VORTAB 

1130 


PHI (3)*PHIT+PHIINC 
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1131 


PHI(9)»PHIT+PHIINC*2.0 


VORTAB 

1132 


R0(3)«RN0T+DELT*2.0 


VORTAB 

1133 


RG(9)«RN0T+DELT 


VORTAB 

1139 


GO TO 53 


VORTAB 

1135 


51 R0(2)»RU(A)-RN0T+.5*DELTAR0 


VORTAB 

1136 


PHI (A ) »A5 .0+.5+PHIT 


VORTAB 

1137 


PHl(2)»-95.0+.5*PHIT 
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1138 


GO TO 53 


VORTAB 

1139 


52 PHIINC»(90.0 + PHIT) /3.C 
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DELT«DELTARL/3.0 
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11 91 


PHI (2 )*PHIINC-9C. 


VORTAB 

1192 


PHI(3)»2.0+PHIINC-90. 
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R0(2)«RN0T+DELT 
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R0(3)*RN0T+DELT*2.C 


VORTAB 

1195 

c 

♦* MAKES THE SECOND HALF OF THE FUSELAGE A 

M TR ROR IMAGE 

V0RTA8 

1196 

c 

+* OF THE FIRST HALF. 
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53 PHI <8 ) =180. O-PHI (2 ) 
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PHI(7)»180.C-PHI(3) 
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PHI (6 ) ■180.0-PHI (A ) 
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R0(8)«R0(2) 
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1151 


R0(7)«R0( 3) 
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RO ( 6 ) “RO ( A ) 
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RETURN 
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1159 


END 


VORTAB 

1155 


C 


VORTAB 

1156 


SUBROUTINE EIRST6(NPAN) 

VORTAB 

1157 

c *+ 

THIS GENERATES THE FIRST 6 CARDS OF THP DATA DECK. 

VORTAB 

1158 


INTEGER WHERE 

VORTAB 

1159 


01 MENS ION T I T L E C 8 ) ,PMACH(7), ALPHA (7) 

VORTAB 

1160 
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C ** A SHORT EXPLANATION OF THE CARD NUMBERS. 

PRINT +,"THE CARD NUMBERS CORRESPOND TO THE CARD NUMBERS" 
PRINT *, "AS GIVEN IN THE DOCUMENT FOR THE VCRLAX PROGRAM." 

C ** CARD NUMBER 1 

71 PRINT *, "ENTER THE DATA FUR CARO NUMBER 1 . •• 

PRINT *, "ENTER TITLE FOR THE AIRCRAFT CONFIGURATION" 

200 READ 100, TITLE 

IF ( EOF (1) >200,201 

201 CALL PRINTC (1, WHERl ) 

IF(WHERE) 42,42,71 

42 WRITE(7,100)TITLE 
100 FORMAT (8A10 ) 

C ** CARD NUMBER 2 

72 PRINT *, "ENTER THE DATA FOR CARD NIIMR.FP 2." 

PRINT *, "METHOD OF SOLVING THE SYSTEM OF FOUATIONS RELATING" 
PRINT *,"BCUNDRY CONDITIONS TO VORTTCITY STRENGTH (ISOLVE)" 

10 PRINT *,"0 FOR GAUSS-SEIDEL" 

PRINT +,»1 FOR PURCELL * S VECTOR ORTHOGONAL I ZATION" 

202 READ (1,*)ISCLVE 
IF ( EOF ( 1 1 1202,203 

203 If (ISOLVE. EQ.C. OR. ISOLVE. EQ.1IG0 TO 11 
GO TO 10 

11 PRINT *,"CH0RDWISE SPACING OF VOR TICES (LAX)" 

12 PRINT *,"0 FOR COSINE LAW SPACING" 

PRINT +,"1 FOR LINEAR SPACING" 

204 READ ( 1, * ) L AX 

IF (EUf ( 1) 1204,205 

205 IF (LAX. EQ.l .OR.LAX.EO.OIGO 10 13 
GO TO 12 

13 PRINT *, "SPANWISE SPACING OF VORTTCFS (LAY)" 

14 PRINT *,"0 FOR COSINE LAW SPACING" 

PRINT *,"1 FOR LINEAR SPACING" 

206 READ (1,* 1LAY 

I F ( EOF ( 1 1 1206,207 

207 IF(LAY.EQ.1.UR.LAY.EQ.0)G0 TO 15 
GO TO 14 
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lb 

REXPAR =0.0 

VORTAB 

1197 


IK IS0LVE.EQ.17G0 70 17 

VORTAB 

1198 


PRINT *, "INPUT THE UV E K -R E L A X AT T ON PARAMETER (REXPAR)" 

VORTAB 

1199 


PRINT *, "INPUT 0.0 FOP MACHINE DETERMINED OPTIMUM" 

VORTAB 

1200 

60 

PRINT *,"0.0 <« REXPAk < .99" 

VORTAB 

1201 

208 

READ ( 1, * ) REXPAR 

VORTAB 

1202 


IF(EUFtl) 7208,209 

VORTAB 

1203 

209 

IF ( REXPAR .GE .0.0 .A .REXPAR .1 T . .99 ) r,n TO 17 

VORTAB 

12 04 


GO TO 60 

VORTAB 

1205 

17 

PRINT *, "HEIGHT ABOVE GROUND (HAG)," 

VORTAB 

1206 


PRINT +, "ENTER 0.0 FOR NO GROUND F F r EC T" 

VORTAB 

1207 

210 

READ (1,*)HAG 

VORTAB 

1208 


IF C EOF ( 1 ) 7210,211 

VORTAB 

1209 

211 

PRINT *, "LONGITUDINAL VORTEX WAKE FLOTATION FACTOR (FLQATX7" 

VORTAB 

1210 

212 

READ (1,+)FLQATX 

VORTAB 

1211 


I F ( EOF ( 1 ) 7212,213 

VORTAB 

1212 

213 

PRINT +, "LATERAL VORTEX WAKE FLOTATION FACTOR (FLOATY)" 

VORTAB 

1213 

219 

READ ( 1, * ) F LOAT Y 

VORTAB 

1214 


I F ( EOF ( 1 7 7214,215 

VORTAB 

1215 

215 

ITRMAX-O 

VORTAB 

1216 


IF( ISULVE.E0.17G0 TO 20 

VORTAB 

1217 


PRINT *, "MAXIMUM NUMBER OF ITERATIONS FOR THE GAUSS-SEIDEL" 

VORTAB 

1218 


PRINT *, "RELAXATION METHOD ( I TR M A X 7 " 

VORTAB 

1219 


PRINT *, " I F 0 IS ENTERED DEFAULT OF 99 WTLL BE USED" 

VORTAB 

1220 

216 

READ ( 1, * 7 I TRMAX 

VORTAB 

1221 


IF C EOF C 1 1 7216,20 

VORTAB 

1222 

20 

CALL PR1NTC(2, WHERE) 

VORTAB 

1223 


IF (WHERE 7 43,43, 72 

VORTAB 

1224 

93 

WRITE (7, 101) ISOLVE , LAX, LAY, REXPAR , HAG, FLOATX, FLOATY, I TRMAX 

VORTAB 

1225 

101 

FORMAT ( IX, I1,9X,I1,9X, I1,BX»4F10.3,7X»T3 ) 
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** i 

CARO NUMBER 3 

VORTAB 

1227 

73 

PRINT *, "ENTER THE DATA FOR CARD NIIMPFP 3." 

VORTAB 

1228 


PRINT *, "NUMBER OF MACH NUMBERS TO R F ANALYZED (NMACH7" 
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50 

PRINT 

VORTAB 

1230 

217 

READ ( 1 , * 7 NMACH 
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IF ( EOF ( 1 7 7217,218 

VORTAB 
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216 IF (NMACH.GT.7)GO TU t? 0 
IF<NMACH.EC.O)GO TC 55 
00 51 J«1,NMACH 

PRINT *, "INPUT MACH NUMBER # ", J 

300 READ <1,*)RMACH(J) 

IF ( EOF ( 1 ) ) 300,51 

51 CONTINUE 

55 CALL PRINTC (3, WHERE) 

IF (WHERE )44,44,73 

A A WRITE ( 7, 103 )NMACH, (RMACHt J ), J«l, NMACH) 

103 FORMAT (I2,8X,?F10.3) 

C ♦* CARD NUMBER A 

74 PRINT *, "ENTER THE DATA FOR CARD MUM R FR A « " 

PRINT *, "NUMBER OF ANGLES OF ATTACK (NALPHA)" 

52 PRINT *,"<■ 7" 

219 KtAD ( 1 , * ) N AL PHA 
IF ( EOF <11 ) 2 19 , 2 20 

220 IF(NALPHA.GT.7)G0 TO 52 
IF(NALPHA.EG.O)GO TU 54 
00 57 J-'l, NALPHA 

PRINT + , "INPUT ANGLE * ",J," IN DFGREFS" 

301 READ ( 1, * ) A LPHA ( J ) 

IF ( EOF ( 1 ) ) 301,57 

57 CONTINUE 

54 CALL PRINTCI4, WHERE) 

IF (WHERE)45,45,74 

45 WRITE (7, 103) NALPHA, ( AL PHA ( J ) , J« 1 , N AL PH A ) 

C ++ CARD NUMBER 5 

75 PRINT *, "ENTER THE DATA FOR CARD NUMBER 5." 

PRINT *, "ASYMMETRIC FLIGHT UR CONFIGURATION FLAG (LATRAL)" 
30 PRINT *,»0 FCR SYMMETRIC FLIGHT AND CONFIGURATION" 

PRINT *,"1 FOR ASYMMETRIC FLIGHT AND/OR CONFIGURATION" 

221 RE AD ( 1 , * ) L A TRAL 
IF ( EOF ( 1 ) >221,222 

222 IF (LAIRAL.EO.O)GO TO 31 
IF (LATRAL .EO. l)GU TU 31 
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GO TO 30 

VORTAB 

1269 

31 

PSI-0.0 

VORTAB 

12 70 


IF (LATRAL .EQ.O>GO TO 32 

VORTAB 

12 71 


PRINT *>" SIDESLIP ANGLt IN DEG&FPS (PST )"• 

VORTAB 

1272 

223 

READ ( 1 > ♦ » P S I 

VORTAB 

12 73 


IF (EOF ( 1) ) 22 3 >32 

VORTAB 

1274 

32 

PRINT *>"PITCH RATE IN D EGR E ES / S FC DNO (P I TCHO ) " 

VORTAB 

12 75 

224 

READ (T> * )PITCHQ 
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1? 76 


IF ( EOF (1 ) >224,225 

VORTAB 
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22 5 

ROLLQ-YAWQ-O.O 
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IF (LATRAL . EQ.OJGO TO 33 

VORTAB 
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PRINT "RGLL RATE IN DEGREES /SFCONO ( ROLL 0 ) " 

VORTAB 

1280 

226 

READ (1,+)RGLL0 

VORTAB 

12 81 


IF(EUF(1) ) 226,227 
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12 82 

227 

PRINT *,"YAW RATE IN DEGREES/ SECOND (YAWOV" 
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1283 

2 28 

READ ( l.,+ ) Y AWQ 
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IF ( EOF ( 1 ) >228,33 
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33 

IF(PITCHQ.NE.O)GO TO 34 
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1286 


IF ( YAWQ . N E .0 ) GO TO 34 
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IF(ROLLQ.NE.O)GO TO 34 
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VINF«G ♦ 0 
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1289 


GO TO 35 

VORTAB 

1290 

34 

PRINT +> "REFERENCE FREE STREAM VELOCITY (VI.NF)" 

V ORTA'S 

12 91 


PRINT *,"IF 0.0 IS ENTERED DEFAULT OF 1.0 IS. USED" 

VORTAB 

1292 

229 

READ ( L, + ) V INF 
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IF (EOF (1) >229,35 
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1294 

3b 

CALL PRINTC (5, WHERE) 
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1295 


IF (WHERE >46, 46, 75 

VORTAB 

1296 

46 

WRITE (7, 104 JLATRAL, PS I, P I TC HQ "L t 0 , Y AWOV VI NF 

VORTAB 

1297 

104 

FORMAT (IX, 11, 8X,5F10.3> 
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** i 

CARD NUMBER 6 

VORTAB 

1299 

76 

PRINT *, "ENTER THE DATA FOR CARD NUM 8 FR ft.." 

VORTAB 

1300 


PRINT *, "REFERENCE AREA FOR FORCF AND MOMENT COEFFICIENTS 

(SREF)" VORTAB 

13 01 

230 

READ ( 1, * > SP6F 
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IF (EOF (1 ) >230,231 

VORTAB. 

1303 

231 

PRINT *, "PITCHING MOMENT COEFFICIENT RFFERENCF LENGTH" 

VORTAB 
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PRINT *, "USUALLY MEAN AERODYNAMIC CHORD LENGTH ( C 8 A R ) " VORTAR 

232 READ (1,*)CBAR VORTAB 

IF (EOF ( 1 ) > 232,233 VORTAB 

233 PRINT *, "X-COOKDINATE UF MOMENT RFFFRFNCF POINT ( X R AR ) " VORTAB 

2 3 A READ (1,*)XBAR VORTAR 

I F ( EOF ( 1 ) ) 2 3 A » 2 35 VORTAB 

235 PRINT ♦ , "Z-CUORDINATE UF MOMENT REFERENCE POINT (ZBAP)" VORTAB 

236 READ <1#*)ZBAR VORTAB 

IF (EOF (1 ) >236,237 VORTAB 

237 PRINT ♦, "TOTAL WING SPAN IN UNITS of SPFF AND C B AR (WSPAN)" VORTAB 

238 READ (1,*)WSPAN VORTAB 

IF (EOF (1 ) ) 2 38 » 2 39 VORTAB 

239 CALL PRINTC (6, WHERE ) VORTAB 

IF ( WHERE ) A7, A7, 76 VORTAB 

A7 WRITE (7, 110 )NPAN,SREF,CBAR,XBAP, ZBAP, WSPAN VORTAB 

110 FORMAT (I2,8X,5F10.3) VORTAB 

RETURN VORTAB 

END VORTAB 

C VORTAB 


SUBROUTINE PR I N TC ( N, WHER E > VORTAB 

C ♦ ♦ ALLOWS A QUESTION TO BE ANSWERED WTTH YF S OR NO, VORTAB 

C ♦♦OR, TRUE UR FALSE, OR ANYTHING LESS THFN 16 CHARACTERS VORTAB 

C ♦♦ THAT STARTS WITH Y , T , N , OR F . VORTAB 

INTEGER WHERE VORTAB 

DIMENSION ANSER(15) VORTAB 

PRINT + , "DO YOU WANT TO RE-EN TFP THF DATA FOR CARD ",N VORTAB 

PRINT ♦, "ENTER Y CR T FOR YES, OP N OP F FOR NO," VORTAB 

A READ (1,1) ANSER VORTAB 

IF (EOF ( 1) ) A, 5 VORTAB 

1 FORMAT ( 15 A1 ) VORTAB 

5 ANS ■ ANS ER ( 1 ) VORTAB 

IF (ANS. EQ.1HTIG0 TC 2 VORTAB 

IF ( ANS. EG . 1HY )G0 TG 2 V0RTA8 

IF (ANS.EQ.IHN)GO TO 3 VORTAB 
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AI’IM NI)T X 15. 


Con 1 itvuod 


If (ANS.I O.IHF )Gll TO 3 VORTAB 

PRINT **"RE-ENIER YOUR ANSWER I TFT JUS T T FT F 0 • " VORTAB 

GO TO A VORTAB 

£ ** YES VORTAB 

2 WHERt * 1 VORTAB 

RE TURN VORTAB 

C ♦* NO VORT AB 

3 WHERE--1 VORTAB 

RETURN VORTAB 

tN0 VORTAB 

c VORTAB 


SUBROUTINE LAST3 VORTAB 

C++ THIS IS USED TO GENERATE THE LAST 3 CARDS «'S 19, 20, £ 21. VORTAB 

DIMENSION X S ( 20 ) VORTAB 

INTEGER WHERE VORTAB 

PRINT**" THE DATA FOR THE LAST 3 CARDS IS INPUT NEXT*" VORTAB 

PRINT *,"THES£ CARDS ARE NUMBERED IP. ?0. AND 21." VORTAB 

C ** CARD NUMBER 19. VORTAB 

PRINT *» "ENTER DATA FOR CARD NUMRFP IP." VORTAB 

1 PRINT + * "NUMBER OF X-S7ATI0NS TCI 0 F F t N E SPATIAL FLOW FIELD" VORTAB 

PRINT *, "SURVEY GRID <21* (NXS)" VORTAB 

bO READ ( 1* * ) NXS VORTAB 

I F ( EOF ( 1 ) ) 5 0* b 1 VORTAB 

bl IF(NXS.GT.20)G0 TO 1 VORTAB 

IF (NXS.GT.O IGG TG 2 VORTAB 

NYS*NZS r O VORTAB 

GO TO 10 VORTAB 

2 PRINT *» "NUMBER OF Y- STATIONS TO D F P INF THF RUTT LINF" VORTAB 

PRINT *# "VALUES OF THE SURVEY GRID (NYS)" VORTAB 

b2 READ ( 1 * * ) N Y S VORTAB 

IF ( EOF ( 1 ) ) b 2 # 5 3 VORTAB 

b 3 PRINT *»"NUMBER OF Z-S1ATI0NS TO DFFTNF THE WATER IINE" VORTAB 

PRINT * * "VALUES OF THE SURVEY GP T D (N7S)" VORTAR 

5A READ ( 1* * ) NZS VORTAB 


13 39 
1 3 AO 
1 3 A 1 
13A2 
13A3 
13 AA 
13A5 
13A6 
I3A7 
13A8 
13A9 


1350 

1351 

1352 

1353 
13 5 A 
13 55 

1356 

1357 

1358 

1359 

1360 

1361 

1362 

1363 
136A 
1365 
13 66 

1367 

1368 
13 69 
1370 
13 71 
13 72 



APPENDIX B. - Continued 



IF<E0F<1> >59,10 
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1373 

10 

CALL PRINTC (19, WHERE ) 
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13 79 


IF(WHLRE) 11,11,1 
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11 

WRITE(9,100)NXS,NYS,NZS 
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100 

F 0 R M A T ( 3 ( 12, 8 X ) ) 
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IF (NXS .GT .0 )G0 TO 9 
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PRINT *> "EXECUTION OF THE PROGRAM T S COM PL E TED" 
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RETURN 
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** 

CARD NUMBER 20. 



VORTAB 

13 81 

9 

PRINT *," ENTER THE DATA FOR CAPO NUMBER 

?0" 


VORTAB 

1382 

5 

PRINT *,"X-STATION VALUES FOR SPATIAL FLOW FIELO 

(XS)" 

VORTAB 

1383 


DO 30 J J ■ 1 > NXS 
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PRINT *, "ENTER X-STATION NUMBER ",JJ 
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1385 

31 

READ <1,*)XS(JJ) 
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1386 


I F ( EOF ( 1 ) >31,30 
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1387 

30 

CONTINUE 
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55 

CALL PRINTC (20, WHERE ) 
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1389 


IF (WHERE >20,20, 5 
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1390 

20 

WRITE (9, 101 ) (XS t J ), j-l,NXS > 
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101 

FORMAT (6F10. 3 ) 
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♦ ♦ 

CARD NUMBER 21. 
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1393 

90 

PRINT *, "ENTER DATA FOR CARD NUMBER >i« 



VORTAB 

1399 


PRINT *, "BEGINNING OF GRID IN BUTT LTNF 

OTRFCTION 

( YNOT ) " 

VORTAB 

1395 

56 

READ (1,+JYNOT 
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I F ( EOF (1) >56,57 
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57 

PRINT *, "Y-SPAC ING OF THF GRID (DFLTAY)" 
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58 

READ (1,+)DELTAY 
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1399 


IF(bOF(l) >58,59 
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Concluded 


59 PRINT "BEGINNING Of- GRID IN WATER LINE DIRECTION (ZNOT)" 

60 READ ( 1# * ) Z NOT 
IKE0M1) )60>61 

61 PRINT *»"Z-SPACING OE THE GR I D ( QEl T A 7 . ) " 

62 READ (1>* )DELTAZ 
I F ( EOT ( 1 ) )62» 63 

63 CALL PRINTC (21>WHERE) 

IF (WHERE ) 41,41* 4C 

A1 WRITE(4> 101 ) YNOT,DELTAY, ZNOT,DEl TA7 

PRINT *> "EXECUTION OF THE PROGRAM TS COMPLETED" 

RETURN 

END 
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16. Abstract 


A method of developing an input data file for use in the aerodynamic 
analysis of axomplete airplane with the VORLAX computer program is described. 
Familiarity with VORLAX input data and panelling methods is assumed. The 
required hardware consists of an interactive graphics terminal equipped with a 
graphics tablet. Software includes graphics routines from the Tektronix PLOT 
10 package as well as the VORTAB program described herein. 

The user must determine the size and location of each of the major panels 
for the aircraft before using the program. Data is entered both from the 
terminal keyboard and the graphics tablet. The size of the resulting data file 
is dependent on the complexity of the model and can vary from ten to several 
hundred card images. After all of the data has been entered, two programs 
READS and PLOTB, are executed which plot the configuration allowing visual 
inspection of the model. 

















